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Chapter 1 
Introduction 



In this work we will discuss decays of the pseudoscalar mesons (tt°, t] and rj') which are 
governed by the chiral anomaly. The chiral anomaly is the non conservation of the axial 
vector current under quantization when gauge fields are present. 

The anomalous decays of the ir° were first measured in the 1960s ([TJ [2] and [3]) and have 
been updated in the recent years [I]. Moreover, many calculations were performed in this 
sector, see e.g. [HI El El IE] and [9]. Some of these decays are of special interest because 
their study permits a deep insight into aspects of modern physics. Through the decays 
7T° — > e + e~7 and 7r° — > e + e~e + e~ the single or double off-shell behavior, respectively, of 
the pion-2-photon vertex can be analyzed. The decay ir° — > e + e~ was first measured by 
[TOT ITT] . Even most recent calculations ([121 ESI EH]) are still three standard deviations 
away from the modern measurements ([IS]). This may open a window for the search of 
new physics. 

By recent measurements of so-called ^-factories, e.g. WASA@CELSIUS, WASA@COSY 
([IS]), CBOMAMI, KLOE@DA<J)NE ([II]) or the CMD-2 Collaboration ([18]), the decays 
of the rj meson have become an important subject of modern hadron physics. Analogous 
to the 7r°-decays, the theoretical calculations of the decays rj — > / + /~7, r] — > and 
7] — > can now be tested by modern measurements. While all the above mentioned 
decays proceed via the triangle anomaly, a study of the box anomaly is possible as well. This 
can be done by analyzing the decays r\ — )■ ti + ix~^ and 77 — > , see e.g. [THJ [20] and 

[2IJ[22], respectively. Particularly these decays admit a study of CP- violations beyond the 
standard model. The decay rj — > 7r + 7r~-f allows a test of a formally unchecked CP-violating 
formalism. In the decay rj — > tt + tt~1 + 1~ CP-violation can be observed via asymmetry 
measurements of the 7r + 7r~- with respect to the Z + Z~ decay planes. Theoretical discussions 
can be found in [23], while measurements were done recently by KLOE@DA$NE (|17j) 
and are in progress by WASA@COSY. 

In the t]' sector experimental data are very scarce and only few theoretical calculations 
were done. New measurements are in progress, but it will take some time until precise 
findings can be expected. Still, theoretical predictions can be made. In addition to the 
analogous decays of the 7/-sector, also the decay 7/ — > 7r + 7r~7r + 7r~ is kinematically allowed. 
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CHAPTER 1. INTRODUCTION 



Since the aim of this work is the study of the anomalous //-decays and analogous decays in 
the 7T°- and ?/-sector, r/'-decays as the above mentioned 7/ — > 7r + 7r~7r + 7r~ process will not 
be discussed here. 

Extrapolated to the chiral point, all anomalous decays are solely determined by the Wess- 
Zumino-Witten Lagrangian (J2U [22]). As we will see for decays into off-shell photon(s), 
we have to take the momentum dependency of the form factors into account. To describe 
that we will apply vector meson dominance (VMD) models. Over the years the so called 
hidden gauge model ([26l[271|28]) has reached good agreement with the experimental data. 
Modern refinements were introduced by the authors of [291 130] and [31] in the last years. 
We will calculate the decay rates and branching ratios for all mentioned decays using the 
hidden gauge model and its modern update in order to discuss the ability of these models 
to describe the experimental data. 

This work is organized as follows. In the second Chapter we will give a short overview over 
pseudoscalar mesons and their decay channels. We will discuss symmetries in general in 
order to analyze the axial anomaly in particular. Thereafter we derive the Wess-Zumino- 
Witten Lagrangian according to [2U and present the contact terms which determine 
the chiral triangle and box anomaly at the chiral point. The concept of vector mesons 
as dynamical gauge bosons of hidden local symmetries will be discussed to explain the 
solutions of the anomaly equation in the presence of vector mesons. Thereby we can present 
the mentioned vector meson dominance models (the hidden gauge model and the modified 
vector meson dominance model [2SJ ED] and [3T]) which will be used in the following 
chapters. 

The third Chapter contains the derivations of the decay rates and branching ratios of all 
mentioned anomalous decays. We will differentiate between the decays proceeding via the 
triangle anomaly and the ones progressing via the box anomaly, because each sector is 
separately governed by closely related form factors. We also present extensive discussions 
of the vector meson form factors for each decay. In the case of the box-anomaly decays, 
we will emphasize the contributions of the CP-violating form factors in addition to the 
leading-order contribution. 

The calculations are presented in detail to be comprehensible for future students. Up to 
Chapter 3, this work can be seen as a handbook for anomalous decays, which is actually 
the reason for the detailed presentation of the calculations. 

In the fourth Chapter we will discuss the results calculated for the decay rates and branch- 
ing ratios with the various vector meson dominance models. We will show the dependence 
of the decay rates on the invariant masses of the outgoing particles and compare our results 
for the branching ratios to other theoretical works and experimental data. 

In the last Chapter we give a conclusion of our work and a brief outlook on further studies. 



Chapter 2 

Theoretical background 



2.1 Pseudoscalar mesons 

In the quark model the different hadrons are classified according to their quark content. 
Because these particles are color-neutral states, hadrons have to be constructed from a 
quark and an antiquark or three valence quarks (or antiquarks). The hadrons constructed 
by two valence quarks , a quark and an anti-quark with color and 'anti-color', respectively, 
are called mesons. The hadrons with three quarks with suitable colors are called baryons. 
These valence quarks give rise to the quantum numbers of the hadrons via their flavor and 
via their symmetry J PC . Here J = L + S is the total angular momentum containing orbital 
angular momentum L and spin S, while P = (— and C = (— 1) L+S stands for parity 
and charge conjugation. Baryons are constructed from three quarks, respectively, three 
antiquarks. Thus they are fermions. Mesons contain a quark-antiquark pair and thus are 
bosons. In the following work we are only interested in light mesons built by up, down or 
strange quarks, which are subject to an approximate U(3) flavor symmetry. The resulting 
nine states can be decomposed into a singlet and an octet state. Written in group notation 
this means: 

3x3 = 8 + 1 (2.1) 
The different mesons can be classified into types according to their spin configurations. 



Type 


S 


L 


p 


J 


J p 


Pseudoscalar meson 













o- 


Axial vector meson 





1 


+ 


1 


1+ 


Vector meson 


1 







1 


1- 


Scalar meson 


1 


1 


+ 





0+ 


Tensor meson 


1 


1 


+ 


2 


2+ 



Table 2.1: Types of mesons 



3 



i 



CHAPTER 2. THEORETICAL BACKGROUND 



S=+l 





0=-i\ 0=«\ fi=+r 



Figure 2.1: nonet of pseudoscalar mesons 



Figure 2.2: nonet of vector mesons 



The different types are shown in Table I2TT1 The nonet of the pseudoscalar mesons (J p = 
0~) and the nonet of the vector mesons are shown in Figure 12.11 and Figure 12. 2[ Here 
the charge increases towards the right and the strangeness increases towards the upward 
direction. Note, that 77 and rj' are not the exact octet and singlet states, respectively. These 
are denoted by r] and 773. The physical, measured particles are mixings between the r] 
and r]s states with an 77-77'- mixing angle 6 mix ~ —20° [32] • These states can be constructed 
from the flavor states according to 



- sin 6> m ix cos 6> mix 
cos miv sin 9 m ix 



(2.2) 



Since we want to study the decay of the pseudoscalar mesons vr°,?7 and 77', we first give 
some general informations of this particles and then show their different decay modes. 



The 7T° is the lightest meson with a mass of m^o = (134.9766 ± 0.0006)MeV [33]. The 
quark content is: 

dd) . (2.3) 



1 



7T° : —= (uu 

Like all mesons, the pion is unstable, however w.r.t. electroweak decays. The decay modes 
are given in Table 12.21 



The 77 and 77' - in comparison to the pion - have also strange quark content: 

r]s '■ —7= (uu + dd — 2ss) , 
v 6 

7/0 : (uu + dd + ss) . 



(2.4) 
(2.5) 



The masses of the mesons are m n = 547.853 ± 0.024 MeV and m v ' = (957.78 ± 0.06) MeV. 
The total decay widths were measured to T v = (1.30 ± 0.07) keV for the 77-meson and 
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Mode 


Branching ratio 


7T° -> 




(98.823 ± 0.034) • 10" 2 


7T° -> 


e^e~7 


(1.174 ±0.035) • 10" 2 


7T°^ 


7 positronium 


(1.82 ±0.29) • 10~ 9 




e + e + e~e~ 


(3.34 ±0.16) • 10" 5 


vr ^ 


e + e~ 


(6.46 ± 0.33) • 10~ 8 




47 


< 2 • 10~ 8 






< 2.7- 10" 7 






< 1.7- 10 -6 






< 1.6 • 10" 6 




V T V T 


< 2.1 • 10~ 6 


vr ^ 




< 6 • 10-4 



Table 2.2: branching ratios of the 7r° decays [M] 



T v i = (0.204 ± 0.015) MeV for the rf. The decay modes and branching ratios are given in 
Table O and Table El 



Mode 



Branching ratio 



77 -> 27 
r] -> 3vr° 
rj — > 7r°27 
rj — > 7r 7r°77 

T] — > 47 

77 — >■ invisible 
77 —7- 7r + 7r _ 7r° 

77 — )• 7T + 7T _ 7 

77 — >■ e + e~7 
7? -> /J + /-i~7 
77 —7- e + e _ 
77 — > 

77 — > e + e~e + e~ 

T) — > TT + TT~e + e~ 

77 — > e + e~fi + fj,~ 
7] — > 

7] —7- /j + /j~7r + 7T 
7] — > tt + it~2j 
7] 7r + 7r _ 7r°7 
77 — > 7r°/_i + /_7 _ 7 



(39.30 ±0.20) • 10~ 2 
(32.56 ± 0.23) • 10~ 2 
(2.7 ±0.5) • 10~ 4 

< 1.2 • 10" 3 

< 2.8 • 10- 4 

< 6 • 10~ 4 
(22.73 ±0.28) • 10~ 2 
(4.60 ±0.16) • 10~ 2 

(7.0 ±0.7) • 10~ 3 
(3.1 ±0.4) • 10" 4 

< 2.7- 10~ 5 
(5.8 ±0.8) • 10" 6 

< 6.9- 10~ 5 

(4.2 - 1.3 + 1.5) • 10- 4 

< 1.6- 10" 4 

< 3.6 • 10" 4 

< 3.6 • 10" 4 

< 2.0- 10~ 3 

< 5 • 10~ 4 

< 3 • 10~ 6 



Table 2.3: branchin; 



g ratios of the 77 decays [M] 



6 



CHAPTER 2. THEORETICAL BACKGROUND 



Mode 


Branching ratio 


rf —> iT +/ K~ri 


(44.6 ± 1.4) • 10~ 2 


rj' — > po7 (including non-resonant ir + ir~^) 


(29.4 ±0.9) • 10~ 2 


rj' —7- ir°ir°rj 


(20.7 ±1.2) • 10~ 2 


rj' — > W7 


(3.02 ±0.31) • 10~ 2 


rj' ->■ 77 


(2.10 ±0.12) • 10~ 2 


7?' -> 3vr° 


(1.61 ±0.23) • 10~ 3 


r/ -> /x + /i _ 7 


(1.03 ±0.26) • 10" 4 


r?' — )• TT + 1T~ U + jl~ 


< 2.3- 10" 4 


rj' — > 7r + 7r~ir° 


(3.7-1.0 + 1.1) • 10" 3 


rj' — > 7r°p° 


< 4 - 10~ 2 


rj' — > 2(ir + ir~) 


< 2.5 • 10~ 4 


rj' — > 7T + 7r _ 27r 


< 2.6- 10~ 3 


r/ 2(7r + 7r~) neutrals 


< 1 • 10~ 2 


??' 2(7T + 7r-)7r° 


< 1.9 • 10~ 3 


?/ 2(7r + 7T-)27r° 


< 1 • 10~ 2 


7/ — )• 3(7r + 7r~) 


< 5 • 10" 4 


rj' —> ir + ir~e + e~ 


(2.5 - 1.0 + 1.3) • 10" 3 


rj' — > e + e _ 7 


< 9 • 10" 4 


7/ —7- 7r°77 


< 8 • 10~ 4 


rj' 47T° 


< 5 • 10" 4 


7/ — )• e + e~ 


< 2.1 • 10~ 7 


7/ invisible 


< 9 • 10~ 4 



Table 2.4: branching ratios of the n' decays [34] 



2.2 Symmetries and anomalies 

Transformations which do not change the physics of a system are symmetry transforma- 
tions. In classical physics this means that the action and thereby the equation of motion 
are unchanged. In a quantum mechanical formulation, e.g. in a path integral formalism, a 
symmetry is given if the Lagrangian and the path integral measure are invariant under the 
respective transformation. The relationship between symmetries and conversation laws is 
expressed via the Noether theorem which says that for every continuous transformation 
that leaves the action invariant there exists a time independent classical charge Q and a 
corresponding conserved current d^J^ = 0. 

There exist many different kinds of symmetries, which are all realized by nature. We give 
a short overview and examples: 

• exact symmetry: examples for exact symmetries are the electromagnetic gauge U(l) 
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or the £77(3) color symmetry of QCD; 

• anomalous symmetry: If a classical symmetry is broken in quantum physics it is called 
anomalous. It is not a true symmetry. An example is the axial U(l) symmetry, which 
will be discussed below; 

• explicitly broken symmetry: the symmetry is explicitly broken by a term in the 
Lagrangian that does not preserve the invariance. If this term is very small it is an 
approximate symmetry. The isospin symmetry is a common example; 

• hidden symmetry: if the Lagrangian is invariant, but not the ground state the sym- 
metry is called hidden. Examples are the broken SU(2)l invariance by Higgs fields 
(spontaneous symmetry breaking) or the SU(2)l x SU(2)r chiral symmetry in the 
strong interactions (dynamical symmetry breaking). 

We are interested in the chiral Ua(1) axial anomaly. The concept of anomalies was in- 
troduced by Adler, Bell and Jackiw ( |35[ 136] ) and also by Fujikawa [37] via path integral 
formalism. We will give a short overview of the calculations given in Chapter 19 of [38] . 
In the massless Dirac Lagrangian the left- and right- handed fermions are decoupled and 
the Lagrangian is therefore invariant under the transformation of the 

$ = e -tfns^ (2.6) 

The corresponding axial current 

hn = (2-7) 

is classically conserved, 

3%, = 0. (2.8) 

This does not hold quantum mechanically when gauge fields are present. The axial vector 
current is built from two fermion fields. Because the product of two local operators can 
induce singularities, we separate their locations x and y, and take the limit (y — x) — > in 
the end. This is visualized in Figure 12.31 

The lowest order contribution (without background gauge fields) results in zero, because 
we have to take the trace over three 7-matrices. The next order contribution instead gives 
a nonvanishing result. Therefore the divergence of the current has the following form, 

0%, = -^e^F^, (2.9) 

which is known as Adler-Bell- Jackiw anomaly [35] 136] . is the electromagnetic field 
strength tensor, F^ u = d^Ay - d u A^. 

Another approach uses the path integral method. The result is that the conservation of 
the axial current clashes with the gauge invariance of the fields. The transformation of the 

1 We use the standard notation of the 7-matrices according to [35]. The parameter 9 is real valued and 
£ nuap j g total antisymmetric tensor in 3+1 dimensions 
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Figure 2.3: higher order radiative corrections of ^>(y)^(x) 



functional measure of the path integral gives an additional contribution via the Jacobian 
J 

D^'D^' = ,r 2 ■ D^D^ (2.10) 
In this case the divergence of the axial vector current reads as 

#M 7 - =(-l) n + 1 r Ml^2-MnP p> (2 11) 

U J5fM \ 1 ) n ! (477)2 r VlH2 ■ ■ ■ r IMn-l^n \* • 1 1 ) 

where n = d/2 with d the number of space time dimensions. Taking d — 4 reproduces 
exactly the Adler-Bell-Jackiw anomaly given in (12 .9|) . The discussion in QCD is quite 
similar. The transformation has the following form 

q = _> q' = e - iTa0 ^ 5 q (2.12) 

where T a = |A a with the Gell-Mann Matrix A a . The corresponding axial current reads: 

The calculations done for the QED sector are valid here as well, either in the point splitting 
formalism via radiative corrections or in the path integral formalizm due to the jacobian. 
This yields the following result 

d^ a = -^- 2 e^G%G% ■ Ti[T a t c t d ] (2.14) 

with the gluon field strength tensor G c and the color matrices t c . Because of the flavor 
trace this term actually vanishes. We see that the axial flavor currents have no anomaly 
in QCD. But there is an anomaly in the electromagnetic sector which is given by 

d^ a = -^e^F^F a/3 ■ Ti[T a Q 2 }. (2.15) 
Here Q is the matrix of the electric charges of the quarks, 

(I 

Q= -\ | (2.16) 
\0 
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Since 

A 3 = [ -1 ] , = -jLe^F^f, (2.17) 

Note, that j^ 3 annihilates a 7r° and we therefore get the anomaly contribution to the decay 

7T° — > 77. 




2.3 Wess-Zumino-Witten Lagrangian (WZW) 

We want to present briefly the effective Wess-Zumino-Witten Lagrangian, which summa- 
rizes and determines the effects of anomalies in current algebra. The discussion will follow 
the presentation of [23] and [25] . 
The QCD Lagrangian is given by: 



m) q 



(2.18) 



with 



C 



d^G u — d v G^ — ig [G At , G h 
(dp - igGJ q 



(2.19) 



where G^ = G^X a /2 is the vector field of the gluons, G^ u is the field strength tensor. 
Because of the inherent nonlinearity and the large effective coupling constant in the non- 
perturbative regime, it is difficult to derive predictions about hadrons directly from the 
QCD Lagrangian. The most successful way to deal with it is lattice QCD, although it needs 
very large computational power and the theoretical insight is limited. In the low-energy 
sector chiral perturbation theory (xPT) as low-energy effective field theory of QCD has 
been successfully applied. Low energy x?T exploits the global SU (3) l x SU (3)r symmetry 
of the QCD Lagrangian in the limit of vanishing quark masses, which was mentioned in 
the previous Section. The lowest order effective chiral action is given by: 



S, 



off 



p 
T 



^xTi[{D^U) (D»U ] )} 



with the chiral unitary matrix 



U = exp I 



(2.20) 



(2.21) 



and f n = 92.4MeV is the physical pion decay constant and P are the pseudoscalar fields 

/>° + 7f* + > - + K+ \ 

K° 



P 



\ 



TV 

K' 



75^ 
K° 



v/3 



'/() 



(2.22) 



Ivs + j= 3 VoJ 
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This Lagrangian does not only contain all symmetries of the QCD Lagrangian, but pos- 
sesses extra symmetries, which are not present in the QCD Lagrangian. According to [21)] 
these symmetries are naive parity conjugation and a U ■<->■ U ' symmetry which counts the 
number of the Goldstone bosons modulo two. These separately are not symmetries of the 
QCD Lagrangian, but form a symmetry (total parity) if they are combined. 
The equation of motion derived from f l2.20p in the case without external fields is given by: 

{firfdpU) = 0. (2.23) 

As shown in [25], the equation of motion which violates naive parity can be constructed by 
adding a symmetry violating extra term with the smallest possible number of derivatives. 
This is given by 

^ (f 2 U^d»U) + Ae^t/ 1 " (d^U) f/ f (d u U) (d a U) U ] (dpU) = 0. (2.24) 

Here A is a constant. Note, that typically there appears a four- dimensional antisymmetric 
tensor due to the violation of the naive parity conjugation. 

The next step is to construct a higher order Lagrangian which leads to this equation of 
motion. The derivation is quite complicated and can be found in [25] . 
Because it is not possible to construct a closed form expression in four dimensions which 
breaks the symmetries the Wess Zumino action is constructed in five dimensions 



T WZ {U) = J d 5 x e^ klm Tr [(($17)17+) {{d 3 U)U^) {(d k U)tf) (($[/)[/*) {(d m U)tf)] 

M 5 

(2.25) 

This integral over a five-dimensional Manifold M 5 can be expressed by Stokes' theorem as 
an integral over the boundary of M 5 , which is ordinary Minkowski space M 4 . 
This action is invariant under global charge rotations U — > U + ie[Q,U], where e is a 
constant and Q the electric charge matrix of quarks. By turning this into a local symmetry 
U — > U + ie(x)[Q, U] also the Wess-Zumino action changes to: 



Q (d u Utf) (d a UU^) (dpUtf) 
f- Q (U^d u U) {U ] d a U) (U^d, 3 U) 



+ 



ie 2 



24tt 2 

x 



J d 4 xe^ (d,A u )A a 

[Q 2 {dpU) U ] + Q 2 U ] {dpU) + QUQU ] {dp) U ] ] (2.26) 
and the effective Lagrangian becomes: 

C = & J d A xTi[(D^U) (iW+)] + nf (2.27) 
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which can be found in many publications, e.g (101 HD H21 S3] and also in 
Note that after expanding U and integrating by parts one can find 



2 

ne 



A = 9^rf° £ ^ F ^ (2 - 28) 

with n = N c the number of colors, which is exactly the part that describes the decay 
7T° — > 77. The author of ref also found that the Noether coupling —eA^J^ describes a 
77r + 7r°7r~-vertex: 

B = - — — UjATr+daTT- V '. (2.29) 

This vertex describes the coupling of a photon to three pseudoscalar mesons and therefore 
the decays 77/77' — > 7r + 7r~7. 

In summary the Wess-Zumino-Witten Lagrangian already determines the triangle anomaly 
sector via A and the box anomaly sector via B. 



2.4 Vector meson dominance (VMD) 

Using the WZW Lagrangian the decay of a pseudoscalar particle into two photons is de- 
scribed very well, as we will see in the next Chapter. Because the theoretical predictions 
for other anomalous decays differ from the experimental data an extended model is needed. 
Vector meson dominance models have reached good agreement in describing the experi- 
mental data of the triangle anomaly sector and are also valid in the box anomaly sector. 
Original work for the so-called hidden gauge model has been done by [26| H3]. Overviews 
can be found in [45] and [46]. We will present this subject according to [27] and discuss 
the so-called 'total vector meson dominance', the 'hidden gauge model' and modern im- 
provements done by [221 ED] and [21] afterwards. 



Vector mesons as dynamical gauge bosons of hidden local symmetries and the 
low-energy theorem 

The nonlinear sigma model based on the manifold G/H = U(3)l x U(3)r/U(3)v is a 
low energy effective theory of massless 3-flavored QCD. The global symmetry G g i b a i = 
U{3)l x U{3)r is spontaneously broken down to the diagonal subgroup Hi oca i = U(3)v- 
The dynamical gauge bosons of the hidden local symmetry H.\ oca i can be modelled as the 
vector mesons (p,u> and <fi). In this framework the KSFR relation and the universally of 
p-coupling can be shown (see [17| 145]). 

Following [27J we introduce the variables £l( x ) and £,r(x) as 

U(x) =£ L (xyU(x) (2-30) 

and non-abelian gauge fields V^(x) = V^(x)T a . The global group U(3)l x U{3)r is gauged 
with external fields and Ar^ such that the corresponding covariant derivatives are 
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given as 

D^ L = (d„ - igV^L + %lA Lii , 

Dptn = (dp-igVjtR + itRARn. (2.31) 
In the following we simplify 

/I 0\ 

A Lfl = A Rfl = eB^Q, Q = -± , (2.32) 

\o o -y 

where is the electromagnetic field. 
Using this, the Lagrangian reads 

£ = La + a ^-V + £gauge fields (2.33) 

with an arbitrary parameter a. In this term £ gaU ge fields stands for the kinetic terms of the 
gauge fields and 



C A = -G-Tr(D^ L -&-D^ R -{;l 
C v = -^Tr [p^ L ■ £t + D^ R ■ 4) 2 • 

By fixing the hidden local U(3) gauge to £r = £' L — e l7T ^ n , the Lagrangians, when expanded, 
yield the following relations: 

vector — photon mixing gy = agf%, (2.34) 
vector — 7r — 7r coulping g V7V7T = ^ag, (2.35) 
vector — meson mass m v = ag 2 f%. (2.36) 

When the parameter a is eliminated from the first two relations the low-energy theorem of 
hidden local symmetry follows, 

— = VI (2-37) 

gv-K-n 

which is known as KSFR relation ([27]). By setting the parameter a = 2 we get another 
KSFR relation namely m v = 2gy nn f% and the universality of the vector meson coupling 
gv-Kn = g- This was postulated by Sakurai in 1960 [19]. 

Solutions to the anomaly equation in the presence of vector mesons 

In the presence of vector mesons the anomaly equation has the form [27] 



2 



5T(Z l ^r,V,A l ,A r ) = -10C-i / Tr 

'm 4 



(dA L f - % ~&A\ ) 



(2.38) 
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with a constants C = —i 2 ^ 2 . Here the gauge variation 5 contains also a hidden local 
symmetry transformation, such that 5 = + $v{ v ) + ^r( € r)- The anomaly equation 

is an inhomogeneous linear differential equation. So the solution will be a special solution 
of the inhomogeneous equation, for which we take the Wess-Zumino action ( I2.26p . and 
general solutions of the homogenous equation. As shown in [271 H5] arid [30] six invariants 
that conserve total parity but violate intrinsic parity can be found. According to [5UJ ETJ 
and [52] two of these invariants are charge conjugation odd and can therefore be omitted. 
So the interesting part of the action is given by the Wess-Zumino action and the remaining 
four invariants, 



wz 



(2.39) 



M 4 



These additional terms do not contribute to anomalous processes like 7r° — > 77 in the 
chiral limit, but can contribute away from the chiral limit via vector mesons decaying 
into photons, since they have the same quantum numbers. So the anomalous Lagrangian 
including all couplings is given by [30J: 

^■anomalous = C-VVP + C-AVP + C-AAP + C-VPPP + C-APPP- (2.40) 

Here V, P and A denote the vector meson, pseudoscalar and electromagnetic fields. The 
different Lagrangian pieces are given by 

N c g 2 



C 



VVP 



C 



AVP 



C 



AAP 



'VPPP 



32tt 2 /. 
N c ge 

N c e 2 

. N c g 



C 



64vr 2 /^ 
. N c e 



C3 ^Tv[d,y v d a V fi P], 
(eg - c^e^ a %A v Tr[{d a V P) Q}P], 
l-cJe^PdpAvdcApTrlCfP], 

c 3 )e^Tr[V,d u Pd a PdpP} : 



ci - c 2 



APPP 



41 - -(ci - c 2 + c 4 )]e^A^Tr[Qd v Pd a Pd^P}. 
24-n A }= 4 



(2.41) 



Here g is the universal vector coupling constant, Q is the quark charge matrix given in 
( I2.32p and Nc is the number of colors. The pseudoscalar fields are defined as 

K+ 
K° 



P = P 8 + P = V2 



7^° 



7T 



V 



7T 



K- 



and the vector meson fields read 



V u = y/2 



( p - 



V2 



K° 



K*+\ 

¥ j 



7^° 



Ivs + j= 3 VoJ 



(2.42) 



(2.43) 
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Note that the prefactors of [30] based on the normalizations of [2H] differ from ours. 
The original terms of the WZW Lagrangian are of course the ones in (I2.4ip Caap and 
C-appp without vector meson couplings. The first three Lagrangians contribute to the 
triangle anomaly, while the last two contribute to the box anomaly sector. 

The various vector meson dominance models 

As one can see, the Lagrangians only depend on three different parameter combinations, 
namely C\ — C2, C3 and c 4 . These parameters have to be fixed by comparison with data. For 
different choices of these parameter sets different vector meson dominance models result. 
Before we present the different models, we can give some general conditions on the pa- 
rameters by discussing briefly the triangle sector. In Ref [30] the amplitude for the decay 
7T° — > 77 is derived from the Lagrangians given in (12. 41ft : 

A(n° ^77) = -i — [ 1_-C4 + C3--3 + ^ ] = -i— [l + 2(c 3 -c 4 )]. (2.44) 

AAP term AVP term VVP term 

To make sure that the vector meson dominance factor is normalized to one at q 2 = and 
to recover the original WZW term we set C3 = c 4 . In this case the second Lagrangian will 
always vanish. We now present the different VMD models. 
The first one is the 'full vector meson dominance': 

c 3 = c 4 = 1 and ci — C2 = -. (2-45) 

o 

Inserting these values of the parameters in (I2.4ip . one can easily see that the second, third 
and last Lagrangian, which are the ones involving photon fields, vanish and so photons can 
only couple to pseudoscalar mesons via a vector meson. Although this model is successfully 
describing data in the triangle sector we will not deal with it in our work for the following 
reasons: (i) In the box anomaly sector at the chiral point, full VMD does not produce the 
WZW interaction. As pointed out in [32], the full VMD model gives a 50% larger value 
at zero four momentum than the anomaly, (ii) The authors of [27] and [1H] calculated 
that the decay rate of the decay u> — > 7r°7r + 7r _ is 2/3 times smaller than the respective 
experimental data. 

The following model, which we present, has the same predictions as the full VMD model 
in the triangle sector, because the values of C3 and c 4 are the same, but it also agrees 
reasonably well with the accepted data referring to the box anomaly as was shown in |32j . 
The parameters are in this case: 

c 3 = c 4 = 1 = ci - c 2 = 1. (2.46) 

Here only the first and last Lagrangian in (12.411) give contributions, while the others vanish. 

We will refer to this model as the hidden gauge model (see e.g. [27]). 

The next model, modified VMD, which is of interest, is a further extention of the hidden 
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gauge model. The authors of Refs. [2§1 ED] and [31] fitted the c-parameters of the different 
Lagrangians and the coupling constant g to the data of the the processes 

— > 7T + 7r~, 

-> (t%)7, 
— >• 7r + 7r _ 7T° 

below a CM energy of 1.05 GeV. Note that the fitted processes do not include rf, so it will 
be of special interest and a nice proof of the extended model to compare the results for the 
different models in the rf sector. 

Ref. [3U] presented various fits, where different data sets have been taken into account. 
Two of these fits were indicated as the ones that represent the data the best. These fits are 
given in table 3 of [30] and were denoted by the conditions (i) Global Fit with ND+CMD 
and (ii) Global Fit with ND+CMD++CMD2. The corresponding values for g, C3 and c± — C2 
are listed in Table [2T4l 





(i) 


(ii) 


9 


5.566 ±0.010 


5.568 ±0.011 


a 


2.364 ±0.011 


2.365 ±0.011 


c 3 


0.927 ±0.010 


0.930 ±0.011 


Cl - c 2 


1.168 ±0.069 


1.210 ±0.043 



Table 2.5: Values of the fitted parameters given by the 'best data sets' in [50] 

Because the fits indeed give very similar values in our calculations we will refer to them as 
one model and state the discrepancies as an error range. 

The vector meson mass is related to the coupling constant g and the parameter a via 
( I2.36p . Therefore there appear different vector meson masses in [30] which vary between 
my = 760 MeV and my = 791 MeV in some of the fits. This will contribute to the error 
of the results in the end. 



e + e" 



e e 
e + e" 
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2.5 Definitions 

2.5.1 The decay momenta 

In the next Chapter we will discuss the kinematics of the different decays of the pseudoscalar 
mesons. Therefore we try to make common definition for the four-momenta. This will 
work for every decay mode except the one into four leptons, which we will relabel in the 
corresponding Chapter. The general definitions are discussed here. Some further relations 
as well as a comparison to the essential kinematics can be found in Appendix A. 
The decay momenta of a pseudoscalar meson P into two, three or four particles Pi are 
defined as follows: 



P(P) -)> Pl (p)+p 2 (k), 

P(P) pf(p + )+p^(p.)+p 3 (k), 

P(P) -> pt{P+)+P2(P-)+Pt( k +)+P±( k -)> 

such that the following relations for the four-momenta are valid in any frame: 



P 
P 



p + k, 

p + + p_ +k, 



P = p + + p_ + k + + fc_ . 

=p = k 

Furthermore, we use the following notations: a four-momentum is denoted by an italic 
letter, e.g. q, whereas the corresponding three-momentum is signalled by a bold-faced 
letter q, i.e. 

( Q°\ 

q \ q x 

q J ' q y 
\ Q Z > 

Moreover, 



q 




(2.47) 



q = q x it + q y y + q z z, 
qy = q z z, 

q ± = q-q|| = q - z (q • z) = q x it + q y y 



(2.48) 
(2.49) 
(2.50) 



where qnis the momentum in the z-direction. 
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2.5.2 The identification of the momenta 

By comparing with the text below the cross section formula (38.19) of Ref. [53] . the fol- 
lowing identifications can easily be made: 

p* = (|p*|, SI*), see below Eq. (38.19) of Ref 
|p* | = |p* ± + P+J = | Pi |, see Eq. (38.20a) of Ref 



Ami 



^s pp m 2 — 2 a/ s pp \l 1 ^ , (2-51) 



i.e., p+ is the three-momentum of particle 1 in the rest frame of particle 1 and 2. Note 
that in analogy to Eq. f)2.5ip 



4777- 

s k k -ml = \^/s^\jl - -y*. = -y/s^ p k . (2.52) 



Furthermore, 



k = (|p 3 |,fi 3 ), see below Eq. (38.19) of Ref [M], (2.53) 
:,|| = | p 3 1, see Eq. (38.20b) of Ref [M] 

_ \ 1/2 (m 2 P , s pp , Skk) _ X 1/2 (s P p,m 2 p,Skk) _ A 1/2 (%>, s kk ,mp) 



2m p 2m p 2m i 



(2.54) 



where 



\(a, b, c) = a 2 + b 2 + c 2 - 2(ab + bc + ca) . (2.55) 



Thus, k is the momentum of particle 3 in the rest frame of the decaying particle. Note 
that Eq. (12.551) allows to relate moduli of momenta in different frames: 

ijjjj _ A 1//2 (?7lp, Spp, Skk) _ \/ s PP ^ l ^ 2 {Spp,m 2 p, Skk) _ \/ s PP gg\ 

2mp mp %\/Spp m P 



The variables s, (3 and A are invariant in all frames. We will give our results in terms of 
these standard variables. 



2.5.3 The identification of the angles 

Let us define the angle 8 P + as the angle between the three-momentum p + and the direction 
k = k|| = z in the P rest frame: 

6 P + = arccos (p+ • k/|p+| ) = Z ( p+, k) . (2-57) 
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Furthermore, we define the angle 9* + as the angle between the three-momentum and 
the direction k = kn = z in the p+p- rest frame: 

9* p+ = arccos • k/|p* |) = Z (p* + , k) = 9 P . (2.58) 

9* + is actullay the angle we will use in our discussions, so we will refer to 9* + as 9 P to make 
things easier. In analogy we can define 0?_ as the angle between the k~ three-momentum 
and the P three-momentum in the k~k + rest frame. Note that the P three- momentum in 
the k~k + rest frame points into the negative ^-direction 

91- = arccos (-kl ■ k/|kl|) = Z (k* , -k) = 9 k . (2.59) 

We also define 0, the azimuthal angle between the plane formed by p+p~ in the P rest 
frame and the corresponding plane formed by k~k + . 



Chapter 3 
Anomalous decays 



In the following Chapter we will discuss the decays of the neutral pseudoscalar mesons 
P G {vr°, 77,77'} that are induced by the chiral anomaly. We differentiate between the ones 
which are governed by the triangle anomaly and the ones resulting from the box anomaly, 
because the structure of the pertinent form factors will be quite similar in the respective 
cases. 




Figure 3.1: triangle anomaly Figure 3.2: box anomaly 



The leading decays induced by the triangle anomaly are discussed next. We add here 
the qualifier 'leading' in order to discriminate these decays from those which involve sub- 
leading sequential decays as, e.g., Bremsstrahlung corrections etc. The discussed decays are 

P -> 77, 
P -> / + /~7, 

p -> i + ri + r, 
p -> i+r, 

where are lepton-antilepton pairs. Obviously only electrons and muons are involved, 
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because the tauons are too heavy. In the case of the box anomaly the tv° decays are dy- 
namically forbidden. The leading decays induced by the box anomaly are: 



P — > TV TV 7, 

P -)■ Tv + Tv-i + r. 



Note that we are not dealing with the decay rj' — > tv + tv~tv + tv~ , although it would be very 
interesting because it was never measured. 

We will present the squared matrix elements and the decay rates of all decays listed above. 
The decays induced by the triangle anomaly will be related to the decay into two photons, 
while the decays induced by the box anomaly are related to the P — > tv + tv~^ decay. We will 
also give explicit expressions for the form factors and the different vector meson dominance 
models which currently are used in the description of the respective decays. 



3.1 P -» 77 

We start with the decay P — > 77. This is probably the most famous anomalous decay, 
because historically it was the first process wherein anomalies were discovered. Expressed 
in terms of the respective momenta it reads: Pp — > 7(e,p)7(e, k), where e and e' are the 
polarizations of the photons. The four momentum of the decaying meson is P = p + k. As 
required by Lorentz invariance, parity conservation and gauge invariance, the amplitude 
has the general structure: 

A(P P -> 7(ei,P)7(e2,A0) = M P (p 2 = 0,k 2 = 0) W e^e^ (3.1) 

The form factor Aip(p 2 = 0, k 2 = 0) holds the information of the decaying particle and 
since the decay products are on-shell photons which are massless it is given by a constant 



M; 



a 



if P = TV 



0. 



^77= cos9 mix - 2V2j^ sin9 mix ^ if P = 77; (3.2) 
a 1 



sin e mix + 2V2j^ cos 6 mix ) if P - 



V ■ 



Here a = e 2 /4Tv « 1/137 is the electromagnetic fine structure constant, f n ~ 92.4 MeV 
is the physical value of the pion-decay constant and f m 1.04/^ and /§ ~ 1.3/^ are the 
singlet and octet Pseudo-Goldstone meson decay constants (see [32]). One can nicely see 
the mixing between 770 and r]s to the pseudoscalar mesons 77 and 77' via the mixing angle 
dmix ~ —20° ([32]). 

Note that the vector meson dominance (VMD) factor is set to unity. In all triangle anomaly 
cases the three different terms Cpaa, £pvv an d £>pva contribute, i.e. the coupling of the 
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pseudoscalar meson directly to the two on-shell photons, the full vector meson coupling 
and the mixed form, where the pseudoscalar meson couples directly to a photon and to a 
vector meson, which decays into a photon. This is shown in Figures 13.31 13.41 and 13.51 






Figure 3.3: direct contribu- 
tion Cpaa 



Figure 3.4: mixed term Figure 3.5: full VMD term 
Cpva C-PVV 



Therefore the amplitude of the r\ — > 77-decay is proportional to 



M P ->■ M P x 




C3 — c 4 

mixed term 



full VMD term/ 



M P x(l + 2(c 3 -c 4 )). (3.3) 



As already discussed the coefficient C3 is equal to C4 in order to recover the usual WZW 
term and to make sure that the VMD term is normalized to 1 in the decay to two on-shell 
photons. So the contribution of the mixed term vanishes. 



3.1.1 Squared matrix element 

The squared matrix element of the decay Pp — > 7(61,^)7(62, k) is given by 

\A(P P j{e,p)j{e, k))\ 2 = \M\\e^ pa E^ v , p , a ,^tfe p 2 k°ei V v ' e{k°' . (3.4) 

If we assume that the polarizations of the photons remain unobserved, the photon polar- 
ization vectors can be summed over. We use the following relation: 

<° tOD = X>(p,A)e M ,(p,A) 

A=l 

= -9w! + - rj , \ 2 - 3.5) 

7] ■ p [7] ■ p) A 

In this formula A labels the two transversal polarizations of the photons, while p is its 
three momentum. Since this expression will be contracted with the four- dimensional an- 
tisymmetric tensor the second and third term cancel out, because they contain either an 
additional p or an additional k. Thus we actually have the simplified form: 

W = -9wf ( 3 - 6 ) 

polarizations 
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This expression can be found in standard quantum field theory books, e.g. equation (5.75) 
of [5B]. Inserting (13. 6p into (13 .4p . one fiends: 

\A(P 77)| 2 = -MpW^pV^^'- (3.7) 

Using p 2 = k 2 = 0, because the photons are on-shell, and the following well-known identity, 
(e.g. equation (A. 30) of [3"8"] ) 

SfMUpa^ U p t a > = -2(gp P '9<T<T' - g P a'9p'a), (3.8) 

we can derive the final expression of the squared amplitude of the decay P — > 77 as: 



|^(P^ 77 )| 2 = M 2 P 2(p ■ k) 2 = Mil (P + k ) 4 

= \M 2 P m%. (3.9) 

3.1.2 Decay rate 

We can now discuss the decay rate of P — > 77. We use the general form for the decay into 
two particles (e.g. (38.17) of [M]) expressed in terms of the squared matrix element. 

dl^ 77 = J--L S \A(P 77)| 2 |ldfi. (3.10) 
2m P lbir z E CM 

Here S = 1/2 is the symmetry factor which appears because of the Bose symmetry of the 
two outgoing photons. Inserting the relations |k*| = E*, because the photon is massless, 
Eqm =m Pi and the squared matrix element (I3.9p . we find the expression of the decay rate: 

Tp^ 77 = -^M 2 P m\. (3.11) 

For the decay ir° — > 77 this expression simplifies to 

2 1 
a m% 

where we inserted Aip = a/(irf n ) from (13. 2p . This is the same as given in [S], namely 

rripgp 

1 7t°— >"n = ~~T7- (o.lo) 

1D7T 

with 

am P M P 

where we just have to use the different form factors given in (13.21) for the respective decay. 
This leads to the results: 
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r 

f?— >77 


■Try'— >77 


theoretical values 


7.73 eV 


0.471 keV 


4.841 keV 


experimental data [M] 


7.7 ±0.6 eV 


0.511 ± 0.003 keV 


4.284 ± 0.245 



Table 3.1: Decay rates with the input data values of [32j(see Eq. 13.21 and the values below) and 
experimental data [31] of the decay P — > 77 

For the decay tt° — > 77 the theoretical value represents the experimental data very nicely. 
It is consistent with most of the results calculated via ChPT, e.g T n o^ iry = 7.78 eV [53] and 
TnO-yyy = 7.74 eV [51]. Other ChPT calculations lead to higher values: T n a^^ = 8.14 eV 

[5B] and Y^a_y ri = [56J. Recent calculations gain the result [57]: ^ifi-^rr = 

7.93 eV ± 1.5%. It is planned by the PrimEx experiment at JLab to reduce the error in 
the experimental value down to 1-2% [58J. 

For the decays t)/t)' — y 7T this is not the case. The reason is the values of [52]. By using 
other values we could be closer to the experimental data. In the following we will give the 
values for the branching ratios with respect to the decay P — > 77 and therefore the terms 
where these values occur would vanish anyway. 



3.2 P -> Z + Z~7 

The next decay that we want to discuss is the one of a pseudoscalar meson P into a photon 
7 and a lepton-antilepton pair , the so-called single off-shell decay or Dalitz decay. It 
is related to the decay into two photons, but in this case one of the photons is off-shell (7*) 
and decays into the lepton pair. Hence the form factors will be very similar, except that 
there is the invariant mass dependence in this case. That is why we will present the final 
result for the decay rate in terms of the double on-shell decay. 

The leptons can be either electrons or muons, but this does not have an effect on the 
kinematical formulae we present. We define the four-momenta for the process Pp — > 
7*(p)7(fc) — > l + (p+)l~ (p-)j(k) so that P = p + k=p + +p- + k holds. Since k 2 = 0, the 
only occuring standard variables are 

dp = 

\(m 2 P ,p 2 ,k 2 ) = \(m 2 P ,p 2 ,0) = (m 2 P -p 2 ) 2 . 




There is only one relevant angle appearing, namely 6 P . This can be taken as the angle 
between the outgoing lepton, l + (p + ), and the pseudoscalar Pp in the l + l~ rest frame, 
namely Q p . 

The amplitude for the decay Pp — > 7*(p)7(fc) —> l + (p+)l~ (p-)'j(k) is given by the following 
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expression: 



A(P ->■ / + r 7 ) = 2M P (p 2 , k 2 = 0)e 



s + ;p-, S-)p u e p k° . 



(3.15) 



Compared to the amplitude of the decay into two photons we see that the polarization of 
the off-shell photon turned into the current j M of the lepton pair, which is given by 

f = eu{k-,sJ)'fv{k + ,s-) (3.16) 

where s± are the helicities of the outgoing leptons l ± . The factor two is a symmetry factor, 
because each of the photons can go off-shell and would therefore turn into the current. The 
polarization of the outgoing photon is e. The form factor A4p(p 2 , k 2 = 0) can be written 
as follows: 

M P (p 2 , k 2 = 0) = M P x VMD(p 2 ). (3.17) 

Here the factor Aip is the one given in the decay to two photons (13. 2p and still includes the 
information of the decaying particle. The factor VMD(p 2 ) is the vector meson dominance 
factor. In this case the diagrams contributing are given by the direct Cpaa term and the 
full VMD term Cpvv, see Figure 13761 and Figure [3~71 As mentioned the mixed term cancels 
because of the factor c 3 — C4 = 0. 





Figure 3.6: direct contribution Cpaa 



Figure 3.7: full VMD term C PVV 



Thus the vector meson dominance factor becomes 

VMD(p 2 ) = 1 - c 3 + c 3 - 



1 



,-r(p 2 ) 

« 2 — 



(3.18) 



Note, that this factor becomes unity in the on-shell case, p 2 = 0, such that the normaliza- 
tion holds. 

In this equations T is the width of the vector meson. It is necessary to consider this width 
in the rj' decays. Otherwise there would be singularities because the 77' mass is larger than 
the vector meson mass. The width is given by 

T(s)=, mv [ — ) |^— j^r) B (8-4*4) (3-19) 
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with g mv = 149.1 MeV, see e.g. [55]. Here we can switch between the different vector meson 
dominance models by inserting different values of C3. In the hidden gauge case (03 = 1) 
the direct term cancels and there is no direct coupling between a pseudoscalar meson and 
a photon. In the case of the modified model (03 = 0.927 or C3 = 0.930, respectively) this 
term will give a small additional contribution. 



3.2.1 Squared matrix element 



In order to calculate the squared amplitude we will use the following projection tensor, 
which is derived in f lA.39j) : 



1/2 1/2 

a_=-l/2 s + =-l/2 

2wy\ (p + -p')ti{p + -p^)^ 



e 2 p 2 x 2 



9m' 



p- 



p- 



(3.20) 



The operator actually becomes Q^ ton given in ( 13. 5 p if the photon goes on-shell. Because 
of the total antisymmetric tensor in the amplitude the second term of ( I3.20p cancels. By 
the same reason only the first term of ( 13. 5 p gave a contribution. The squared amplitude 
then reads: 



Ae 2 M 2 P \VMD(p 2 ^ 2 
Ae 2 M 2 P \VMD(p 2 ^ 2 



\a(p / + /- 7 )r = - ^;: 2 ;r v " 71 w^o ff, (p-,p + )o&oy fc '/ fc8 ' 



(p 



(p 



2\2 



" fjbu pa^ pi V p' a' 



I'll 



, (p + —p Y(p + —p y' 



p 



-g uu ')p p k a p p 'k a ' . 

(3.21) 



The structure gets simpler when the following momentum relations are inserted and the 
antisymmetric structure of the epsilon tensor is used again: 



k = P — p 
The squared amplitude is then: 



P-=P~P+ 



(3.22) 



\A(p ^/ + r 7 )| 2 

4e 2 M 2 P \VMD{p 2 )\ 2 

(p2) 2 



(3.23) 



We can now switch to the rest frame of the pseudoscalar meson where the relation P^ = 
mp5 M ° holds. Note the sign change due to g n = —5 n . Thus the squared amplitude reads 
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now: 



4 e 2 M 2 P \VMD(p 2 ^ 2 
2e 2 M 2 P \VMD(p 2 )\ 2 



(p 



,2^2 



m. 



£ ijk £ ijk> p k p k> _ tij* e i?Vih*W 
2 P dp+ P sin p 



(3.24) 



The squared matrix element of the decay P —t Pi 7 can now be given in terms of the 
standard variables as 



\ A (P ;+r T )| 2 = Ae2M2plVMD{Spp)l2 mUm 2 P - s pp ) 2 [2 - (3 2 k sm 2 9 k ] . (3.25) 
3.2.2 Decay rate 

In order to calculate the decay rate we have to deal with the phase space first. The phase 
space for a three-body decay can be found, e.g., in (38.19) of [3lj : 

d$ 3 (P P ;p + ,p_,/c) = ^ ^ | P ;||k|dm p+p _dfl; + df2 fc . (3.26) 

In our case this leads to the following form of the decay rate: 

dr = 7^31 32m \A\ 2 </s^p~/3kEy d^Sp^ dcosflfe dcos6> p d0 d0 7 . (3.27) 

After integration over the angles we find the following result 

dr-^fl=f^K-*)V,[.-fl^ (3,s) 



where we inserted the relations 



2 2E 7 

Sp„ = TTlp(l ), 

mp 

l^&J~Sp~p~ = hds P p = —nip&E^. (3.29) 



J PP^V PP 2 PP ""^ uiJ 7 

Note that the first factor is exactly Tp^ 77 = ■ The final expression for the decay 

rate of the decay P — > / + /~7 can be expressed in terms of the decay P — > 77 as 

dT(P -+ J+r 7 ) = Tp^ al \ MD 6 {Spp)l2 (m 2 P - s pp ) 3 p [3 - /3 2 ] d, pp . (3.30) 
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3.3 p -> i + ri + r 

We will now discuss the double off-shell decay of an pseudoscalar particle P in two lepton 
pairs. The decay is again related to the P — > 77 decay, but in this case both of the photons 
go off-shell (7*7*). The general structure of the form factors will be very similar, but there 
will be a dependence on the invariant masses of both lepton pairs. 

The leptons can be either muons or electrons, so that we are dealing with the following 
three decays: 

(i) P ->■ fi + fi~e + e~, 

(ii) P — > ji + jjT /i + fjT , 
(Hi) P — > e + e~e + e~. 

In the case that the pseudoscalar meson is a 7r°, only the decay into two electron pairs (Hi) 
is possible, of course. 

We define the momenta of the decay P(P) —> If (pi)li (p2)l£ (Ps)^ (P*) as follows: 

P = Pi + P2 + PZ + Pi , 
Pij = Pi +Pj- 

The relevant variables can now be written as: 



Sij = (Pi+Pjf, 




with rriij = mi is the lepton mass of the outgoing leptons i and j. We define the angle 
Qij as the angle between and pp in the l(pi)l(pj) rest frame, the angle as the angles 
between the planes formed by l(pi)l(p2) and l(pz)l(p±) and as the angles between the 
planes formed by l(pi)l(p±) and l(p2)l(pz)- 

3.3.1 Amplitudes 

Calculating the decay of a pseudoscalar meson P into two lepton pairs we have to differen- 
tiate between the decay into two different lepton pairs (i) and the decay into two identical 
lepton pairs (ii,iii). The matrix element squared |Al| 2 of the decay (i) can be calculated 
relatively straight forward by using the invariant decay amplitude (see Figure [378]) : 

Ai(P fi + fi~e + e-) = ^^e^ upa ( Pl + p 2 ) u (p 3 + PaTHp2)Yv( Pi ) -^(p^Yiips). (3.31) 

S 12 S 34 
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Figure 3.8: Double Dalitz diagram for A\ Figure 3.9: Double Dalitz diagram for A 2 

For the decays into two identical particle pairs (ii, in) we have to take another amplitude 
A 2 into account (see Figure l3~9l) : 

a 2 {p -> i + n + r) = -^Ple^ pa { Pl +p A )%p 2 + p 3 y^p 4 )Yv(pi)-u(P2h p v{p 3 ). (3.32) 

S 14 S 23 

In this case the a total amplitude is A = A\ + A 2 . In order to give the squared matrix 
element \A\ 2 it is necessary to calculate the direct term |v4i| 2 , but also the crossed term 
\A\\ and the interference term IRei^AxAty. The crossed term will look the same as the 
direct term when the variables are exchanged. To determine the decay rate we will have 
to integrate over the variables. Therefore the direct and crossed term will give the same 
contribution to the decay rate. 



3.3.2 The reaction P — > fi + fi e + e 

The invariant decay amplitude mentioned above can also be written in the following form: 

Ai(P -> /i + /i"e + e") = -^^e^j^Jp 3 ,p A )p^J^ ) (pi,p 2 )Pu . (3.33) 

P12P34 

Here j? ef * and j?^ are currents of the respective lepton pairs. In order to calculate the 
squared amplitude we need the projection tensor, which we used in (I3.20p . see also f lA.39j) . 
Here it reads 



9W „2 } „2 



0^{pi,p 2 ) = e 2 p\ 2 x 2 

P12 J P12 



(3.34) 



The squared matrix amplitude then reads: 



\A^{P /i + /i-e + e-) = J^e^V-'^'^ , (pi J p 2 )^ 2 0^(p 3 ,P4)j4^ (3-35) 



Si oS 



12°34 
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The following steps are as usual. We insert the relations between the momenta to make 
use of the antisymmetric tensor: 

p 2 = P-pi-pu, 
p 3 = P-P4—P12, 
P34 = P-'Pi2- 

We can now switch to the rest frame of the pseudoscalar meson, where we use P^ = mpS^ 
and g %% = —5 n . The squared amplitude can now be expressed in standard variables as 
follows: 



e 4 |M| 2 



S12S34 
e 4 \M\ 2 
S12S34 



A(mp, S12, s 34 ) 
A(mp,si2,s 3 4) 



2 - (3 2 l2 sin 2 12 - fill sin 2 # 34 + /? 2 2/ 3f 4 sin 2 12 sin 2 # 34 sin 2 
1 + (1 - (3 2 12 sin 2 12 )(1 - /3 3 2 4 sin 2 # 34 ) J sin 2 



2 — /3 2 2 sin 2 ^12 — /?34 sin 2 # 34 ] cos 



(3.36) 



The result is symmetric under the exchange of the momenta p%, p 2 and p%, p 4 . 

Note that the last line agrees with Eq.(16) of [8], expressed in standard variables. It also 

agrees with the result given in Appendix B of [7]. 



3.3.3 The reaction P -> l + l~l + l~ 
Direct and Mixed Term 

We will now calculate the squared matrix element of the decay into two identical lepton 
pairs The part of |Al| 2 is of course the one given in f )3.36p . We can get \A 2 \ 2 by 

replacing p 2 and p 4 : 



\a 2 \ 2 (p ->• i+ri+i- 



S14S23 
e A \M\ 2 
S14S23 



X(m 2 P , s u , s 23 ) 
X(m 2 P ,s u ,s 23 ) 



2 - /3 2 4 sin 2 6 M - /3 2 3 sin 2 6 23 + /3 2 4 /3 23 sin 2 14 sin 2 23 sin 2 
1 + (1 - (3 2 lA sin 2 14 )(1 - (3 2 23 sin 2 23 ) ) sin 2 



+ 2 - /3£ 4 sin^ 14 - /3 23 siir 23 cos 



(3.37) 



Note that this formula is given in different variables than used in all other decays. To 
compare this terms it is necessary to express the new variables in terms of the ones used 
before. The translation formulae are given in Appendix B. 
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Interference Term 



We will now present the term 2Re(*4.i*4. 2 ). The calculation is much more extensive than 
the calculations presented before, so we will give less steps than before. 
We have to deal with the traces at first. In the calculation there are traces over products 
of four, six and eight 7-matrices. We therefore use the formula 



Tr(^ 2 ...^ 2n ) = a x ■ a 2 Tr(fi 3 ..4 2n ) - a x • a 3 Tr(^ 2 ...^ 2n ) + ... + ai • a 2r JV(;i r .^ 2re _ 1 ), 



which can be looked up in standard books on quantum field theory e.g. |60j . 

This leads to very long calculations in which errors, especially in signs, can easily occur. 

Hence the following result was rechecked with the algebraic program FORM (see [6Tj). 



^\M\ 2 e lu/piT e ll i v ipf a >g pp 

2222 
P12P34P14P23 



"(pi +P2Y(P3 +P4T(P1 +PaY'{V2 +P3T' 



X 



m 2 + pi ■ p 2 ) (p£p% - p£ p%j + (m 2 + p 3 ■ p 4 ) \p$p$ - p^ p% 
(m 2 - p 2 ■ p^ (pip'i + P1P3) + O 2 - Pi ■ Pa) (P2P4' - P2P4 
(m 2 + p 2 ■ p 3 ) (piPa - P1P4) + (m 2 + Pi • Pa) (p 2 P 3 ~ P 2 P3 



(3.3* 



This equation is invariant under the permutation of p 2 and p^. As one can easily check, 
the first and second term will turn into the sixth and fifth term, while the third and fourth 
term stay invariant. 

Note that this is not the formula given in Appendix B of [7], which has an typological 
error. This was recognized by the authors of Ref . [8] , too. 

Next we have to deal with the four- dimensional e-tensors. Here we have to contract the 
e-tensors explicitly. We will again use FORM to control the calculation. Inserting the 
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relations P2 = P12 — Pi and p 3 = p 34 — Pa and identifying p 2 = m 2 we find: 
e 4 |M| 2 



\A ■ 



x 16m | (p 12 ■ p 34 ) - p 12 p u 
+4m 2 • pi [(pis • p 34 ) 2 - P 2 2 /4] - 2 p 2 2 P 34 Pi 2 ■ P34 

- 2^12 • PiPl ■ P3A [p\ 2 + pIa + 4j?12 • P34] + 2(p 12 ■ P4) 2 P34 + 2 G°1 ' P34) 2 p? 2 

+ P12 • Pi [(Pn - Pu)Pi2 ■ P34 + ^PnPu] + Pu ■ Pi [(p 34 ~ Pu)Pi2 ■ P34 + 4p 2 2 p 34 ] 



+2< -8pi ■ p 4 Pl2 ■ PiPl ■ P34P12 ■ P34 - 2^12 • P4Pl ■ P34 [PnPu ~ 2 (Pl2 ' P4) 2 ~ 2(pi • p 34 ) 2 ] 



+ 4p! ■ p 4 [2p 2 2 ^ 4 (p 12 • p 4 + Pi ' P34) - Pl2 ■ P34 ((Pl2 • P4) 2 + (Pi • P34) 2 ) - P 2 2 P34(Pl2 ■ P34 

1 

2 



8(P1 • P 4 ) 2 [(Pl2 • P34) 2 - P12P3J - ^12^34 [(Pl2 • P4) 2 + (Pi • P34) 2 ] + 7^2^34 f 1 • ( 3 - 39 ) 



All terms of this formula are still valid for any Lorentz frame. So we can make use of the 
invariants given in Appendix A. 2. This leads to another long calculation, but a check by 
any algebraic program (e.g. Maple) shows that this formula, written in standard variables, 
is given by 



e'\M\ 2 X(m%p 2 2 ,p 2 u ) 



|2 



P12P34P14P23 16 
x \ 4si2/3i 2 2 s 3 4/334 2 sin 2 12 sin 2 34 cos 2 4> + 64m 4 

s! 2 "/3i2v / ^34/ 3 34 sin 12 sin 34 cos 1 



L4 



32m 2 - (m 2 - p 2 2 - p 2 4 ) (/3 12 cos 9 12 - /3 34 cos 6> 34 ) 
- 2p 2 2 p 34 (l - /3 12 f3 34 cos 6>i2 cos 34 ) (2 - /3 2 2 cos i2 - /3 2 4 cos 34 ) 

+ 8m 2 (/3i2 cos 6»i2 - f3 34 cos 34 ) (p 2 2 /3i 2 cos 9 12 - p\ 4 $34 cos 34 ) j . (3.40) 

Note that this formula is exactly the one given in [8] . Only the definitions of the angles 
differ, as mentioned in the Section referring to the kinematics (see Appendix A). 
The s 23 and S14 terms in the denominator can be expressed again via the relations given 
in f lRTj) and (TB^jl . 
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3.3.4 Decay rate 



The decay rate can now be stated in terms of the four-body phase space which is given in 
(IA.46P in a general form: 



dr = \A\ 



P12 A34A 1/2 0i2, m 2 P ,s u ) 



m% ■ 2 



15 



7T 



dsi2 ds34 dcos 9 12 dip dcos #34. 



(3.41) 



Inserting the squared matrix element (13. 36p and integrating over the angles, we get the 
final expression for the decay rate of the process P — > fi + fi~e + e~: 



dTi 



e 4 \M(s 12 , s u )\ 2 (3 12 f3 34 \2 (m 2 P , s 12 , s u ) 



nil 



2 io . . Sl2 . S34 



1 " ~ (/5i2 + &) + \ft2P, 



dsi2 ds 34 . 
(3.42) 



In all decays induced by the triangle anomaly we will later give the branching ratios with 
respect to the decay rate of P — > 77. This simplifies the calculations. Moreover, they 
will also be more precise, because many terms of the form factors are the same. Thus we 
neither have to take care of the various values of the rj /rj s mixing angle nor the coupling 
constants fo and /§. In this case the branching ratio simplifies to: 



dTi 

dr 



77 



2a 2 /3i 2 /3 34 A 2 (m 2 P , s u , s u ) 
9tt 2 ■ S12 • s 34 
1 



dsi 2 ds 34 x |VMDi(si2,S34) 



(3.43) 



Here the factor VMDi represents the vector meson dominance model that we used. 
The interference term in this case becomes: 

£12 = a 2 \(m 2 P ,p 2 2 ,p 2 M ) 



77 



4vr 3 S12S34S14S23 

^12 S 34/3 3 4 '• 



4si2/3i 2 s 34 /3 34 sin 2 9 12 sin 2 6*34 cos 2 



x Re(^MDi 2 (si2, s 34 , s u , S23)) 
64m 4 

32m 2 — (m 2 



- v / ^i2/ 3 i2v / ^34/ 3 34 sin 12 sin 6> 34 cos 

- 2p 2 2 p 2 4 (l - (3 12 (3 34 cos 9 12 cos 34 ) (2 - /3 2 2 cos 9 12 - /3 2 4 cos 6> 34 ) 

+ 8m 2 (/3i2 COS 6>i2 - /3 34 COS 34 ) (pl 2 Pl2 COS 012 - P 34 /?34 COS 34 ) | 

In this case the vector meson dominance factor depends on all four variables. 



Pl2 ~ Pu) (Pl2 COS 012 - /?34 COS 34 ) ' 



(3.44) 



3.3.5 Vector meson dominance factor 

In the final step we will show the terms of the used vector meson dominance model. We 
have again contributions from the terms Cpaa and Cpvv, while the £pvA-term vanishes, 
because of the prefactor C3 — c 4 . This is shown in Figure 15-lUl and Figure 13". Ill 
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i+ 



Figure 3.10: direct contribution Cpaa 



Figure 3.11: full VMD term C PVV 



Thus the factor for the direct and crossed term can be constructed as: 



VMD^s^su) = 1 -C3 + C3- 



mi 



mi 



S12 - im v T(si 2 ) m v - s 34 - im v r(s 34 ) 



(3.45) 



The normalization to unity can be tested by setting the invariant masses to zero. 

In this equations, T is again the width of the vector meson as given in f)3.19p . 

The corresponding factor for the interference term looks somewhat different. We build it 

up via the expressions VMD(su, S34) and VMD*(s\4, s 23 ). Thus we have to calculate the 

real part of the term 



VMD 12 (s 12 ,s M , s u , s 23 ) 



1 - c 3 + c 3 - 



mi 



in 



S12 - im v T(s 12 ) my - s 34 - im v r(s 34 ) 



m 



v 



mi 



mi 



x 



m v ~ s 23 + im v r(s 2 3) m v - s u + im v T(s 1A ) 



(3.46) 



Note the dependence on the four variables S12, s 3 4, S14 and s 23 . The variables S14 and s 2 3 
can be expressed by the usual variables as reported in (IB.ip and (IB.2p . 
In the discussion of the results we will compare the values of two different vector meson 
dominance models, the 'hidden gauge' and the 'modified' one, with the values calculated 
without a VMD-factor. The different models can be specified by setting c 3 = 1 (hidden 
gauge) and c 3 = 0.927, respectively, c 3 = 0.930 (modified VMD model). The two values 
for c 3 of the modified model will give approximately the same results, so we include them 
via the errors. 
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3.4 P -> /+/- 

The final missing decay that proceeds via the triangle anomaly is the one of P(q) — > 
l + (p')l~ (p) ■ A lot of work was done on this subject already by several groups e.g. [Ml E31 
l64"t l65| 166] or [67]. We will mainly follow the work of [T2| 168] . There the dependence on 
many different form factors was discussed, so we can study how the new optimized vector 
meson dominance model acts in comparison. In [13] additional corrections were handled 
and the results were completed. 




Figure 3.12: The amplitude of P ->• l + l~ 



We define the momenta of the decay P(q) — > l + (p')l (p) such that the following relation 
holds: 

q=p + p'. (3.47) 



We use the variable fii{q 2 ) — dl — and of course the momenta and masses of the 
respective leptons to describe the decays. 
According to [67J the Feynman amplitude is given by: 



e 1 



M{P^l + l ) = j- u(p, sJ)Av{p\ s. 



e 1 



where 



Tr [v(p',s+)u(p,S-)A] (3. 

Jtt 

Here the form factor VMD(k 2 , (q — k) 2 ) contains the vector meson dominance input and 
will be discussed later in detail. 

Using the informations given in Appendix C we can write the branching ratio of the decay 
p — > l + l~ in terms of the reduced, dimensionless amplitude A: 

/ \ 2 

P^l+l- r, I am l \ o rJ2\\ At 2 \ i2 



r 

1 77 



/ ami V 
\ivmp J 



^(q 2 )\A(m 2 P )\ 2 (3.50) 
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with 

(3.51) 

as given in [55] and [T2]. Note that [57] differs from other calculation, e.g. [55] by a factor 
2 in the form factor. 

To solve this integral we follow the dispersion approach applied in many publications before, 
e.g. [621 E21 ESI [12]. Therefore, we will first calculate the imaginary part of the amplitude 
A using the Cutkosky rules [69], [70] (for further calculations see Appendix C): 

^ - m^m- (3 - 52) 

This is the contribution of two on-shell photons in the intermediate state. Therefore, the 
form factor is trivial and model independent, i.e. VMD(0,0) = 1. We can give a lower 
limit for the reduced branching ratio by using only the imaginary part (|*A| 2 > IuuA 2 ). 
These unitary bounds are given in Chapter 4. 

In order to calculate the branching ratio we have to deal also with the real part of A. 
According to [65] the once-subtracted dispersion relation is given by: 

Re^(g 2 ) = A(q 2 = 0) + ^ f° ds . (3.53) 

Jo s{s - q 2 ) 

The integral can be calculated using analytical programs. We are only interested in the 
integral for q 2 > Am 2 Q, where mi is again the lepton mass, which is given by the following 
expression [12] : 



Ro4(g 2 ) = A{q 2 = 0) 



A(g 2 



i + A(g 2 )7 12 V i + 



(3.54) 



with the dilogarithm function Li 2 (z) = — J Q z (dt/t) ln(l — t). To the leading order in 
(mi/mp) 2 this is given by 

ReA(m 2 P ) = A(q 2 = 0) + In 2 ( — \ + (3.55) 

\mp J 12 

Thus all the nontrivial dynamics of the process are contained in the subtraction constant 
A(q 2 = 0). The derivation of this constant is summarized in Appendix C according to the 
procedure of Ref. [12]. To the first order in m; /my, where my is the vector meson mass, 
which appears in the calculations given in Appendix C, the expansion is given as follows!: 

A{q 2 = 0) = 31n + X M (3.56) 



1 The expressions for q 2 < and < q 2 < Amf are given in Rcf. 
2 Higher order corrections can be found in Ref. [14] , 
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where \i is a scale parameter, which can be set to my. The constant x{p) depends on the 
VMD factor in symmetrized Euclidean kinematics: 

5 3 f°° j i ( 1 \ dVMD(-t, -t) 



5 3 f°° ( t \ 

xM = -i + 2j dHn UJ 



, it VMD(-t,-t)-l + it VMD(-t,-t) 



■2 



(3.57) 



The VMD factor now depends only on one invariant mass t which simplifies the calculations. 
Note that A(q 2 = 0) is independent on the parameter fi after integration of t. Also note 
that VMD is expressed in Euclidean space-time and we therefore do not need a width, 
because there are no singularities. The form factor can be written as follows: 



VMD(—t, -t) = l-c 3 + c 3 2 V 2 V (3.58) 

m v + t m v + t 



where c 3 is again the c-parameter that differentiates between the various VMD models. 
In the case of 7/ decays there is an additional contribution to the imaginary part ( [14] ) . 
Because the mass of the 7/ is greater than the ones of the u and p mesons, one vector 
meson can be on-shell. This contribution can be calculated again with the help of the 
Cutkosky rules. The additional contribution is given as: 

7T / 1\ 2 , /l-P 



AbnA = -p V --) iB^jOfz-l) (3.59) 

where z = {m p /my) 2 . Note that in comparison to [H] the numerator and the denominator 
are exchanged. 



3.5 P 7T+7T 7 and P ->• n + 7r l + l 

In the next Section we will deal with the decays P — > 7r + 7r _ 7, P — > n + 'K~e + e~ and 
P — > 7r + 7r _ /.i + /i _ . The three decays are very similar in their basic structure, since the 
photon in the P — > 7r + 7r~7 decay can be replaced by an off-shell one that decays into a 
lepton pair. Because of kinematic reasons only the 77 and rj' decays are possible, but not 
the 7T° decay. All of the decays governed by box anomaly proceed as shown in Figure 13.131 

These reactions are quite interesting, because there is the possibility to measure CP viola- 
tion as in the equivalent kaonic decays [7T] . To model this we have to take the usual mag- 
netic vertex with a form factor M into account and add another vertex, which is constructed 
as an CP violating electric dipole operator. The corresponding form factor E was intro- 
duced by [72] for the P — > 7r + 7r _ 7 decay and modified by [23] for the P — > 7r + 7r~/ + /~ decay. 
The form factors will be presented after we have handled the kinematics of the respective 
decays. The structure of the decay amplitude is visualized in Figure 13.141 Note that the 
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Figure 3.13: box anomaly 



contribution of the Bremsstrahlung term is negligibly small (BR(rj — > tttt) < 3.5 • 10~ 14 ), 
see [73], as otherwise the electric dipole moment of the neutron would have been measured 
with a value bigger than the experimental upper bound tabulated by PDG 




7T _ 

M\ (CP conserving) 



P 





£i (CP violating) 



Figure 3.14: CP conserving magnetic (Mi), CP violating electric (Ex) and Bremsstrahlung 
contribution 



In the squared amplitude there appear terms proportional to |M| 2 and \E\ 2 which are 
both CP conserving. The relevant terms for a CP violating contribution are proportional 
to ~Re(ME*). As we will see it is not easy to construct these terms and even harder to 
measure them. 

We will use the definitions of Ref. [72] for the four-momenta of the decay P(P) — > 
n + (p + ) n~(p_) 'j(k), such that the following relation is valid in any frame: 



P = p + + p_ +k , 



(3.60) 



whereas the four-momenta for the decay P(P) tt + G»+) tt (p_) l + (k + ) I (&_) are defined 
according to Ref. [23] as 

P = p++p + k + + k. . (3.61) 



=k 



In order to make the kinematics more obvious we will rename the standard variables s 
(3 p and 9 p as s nn , (3 n and 9 W labelled by the respective particles and not by the momenta. 
Of course s^k turns into s ee and so on for the P — > 7r + 7r~/ + / - case. 
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3.5.1 The reaction P — > ir + ir 7 

We start out with the decay P — > n + ii~ r y and therefore with the invariant decay amplitude 
given in Eq. (2) of Ref. [72]: 

A (P _> vr + vr ^ 7 ) = _L (M G e, va ^k v p%p p _ + £ G [(e ■ p+)(k ■ p_) - (e ■ p_)(A; • p+)]) . 

(3.62) 

Here we can nicely see the usual magnetic form factor M attached to the total antisym- 
metric tensor, which is typical for anomalous decays, and the new contribution of the CP 
violating electric form factor E. Note that the form factors of Ref. [72] and Ref. [23] differ 
in structure and normalization. Therefore, the subindex G is added here. 
The following steps are straight forward. We first insert p_ = P — p + — k, so that we can 
make use of the antisymmetric e-tensor. The amplitude then reads: 



A ( P ^ 7r + 7r - 7 ) = _ M G e, ua ^k»p a + pP + E G [(e • p + )(k -P)-{e-P){k- p + )\ 

VTLp \ 

(3.63) 

In the P rest frame, where P M = mp5 M °, the amplitude is then given by 



A(P^ 



7T + 7T 7) 



m. 



ml 



E 



G 



(-e-p + )(E 1 



M G k- (eAp+) -E G e-p. 



(3.64) 



The squared amplitude reads 



|-4| 2 (P -> 7T + 7r"7) 



E: 



ml 



M G | 2 |k- (eAp +± )| 2 + \E G \ 2 \e ■ 



+E* G M G k-(p +± Ae) (e • p+jj* + M G E G k • (p +J _ A e 



l e • P+±, 



In Ref. 



the following polarization vectors are used: 

(p+ A k) A k k(p + ■ k) - p_ 



£2 



|(p+Ak)Ak| |k(p+-k)-p 
p + A k p +J _ A k 



|p+Ak| |p +± Ak| 
Thus the unpolarized squared decay amplitude reads 



kA(-p +± ). 



J2\A\ 2 (P^^ 

pol=l 



Ti 



7) 




(3.65) 
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The contributions of the mixed term vanished when both photon polarizations had been 
summed. This means that the CP violation in this decay cannot be found, if the polariza- 
tions of the photons are not measured explicitly. 

The squared amplitude f !3.65j> . expressed in standard variables s nn and 9 n , reads 



E, A i2/ a \ H m pi s 7T7r, O)^/? 2 sin 2 6 n f 2 2 
|^P^7r+7r- 7 | (^7™?^ = — g lM W G| + \&G\ 



pol=l 



In terms of J5 7 and 



(3.66) 



4m! / 4ml 



s TT (£' 7 ) y mp — 2mpE, 



n = 

(see Eq. ( I3.69P ) it is given by 

2 _ £ 2 (l-2£ 7 /m,)/3 2 sin 2 ^ / x 
J]|-4p^r +7 r- 7 | 2 (£ 7 A) = ^ j^ 1 ( \M G \ 2 + \E G \ 2 ). (3.67) 

pol=l mp ^ ' 

3.5.2 Decay rate 

The decay rate for the three-body decay P — > 7r + 7r~7 is given by relation (38.19) of Ref. 

1 1 

(2^)5 16m : 



= Tc~ 2 l^| 2 Ip+*I 1^1 drrimr dcos^jr dcj)* dcos# 7 d0 7 



/ s. 3 3/2 

1 - % ^^sin 2 ^ {\M G (s^ } \ 2 + \E G \ 2 ) dy^dcos^ 
v Trip j Trip 



2 12 7T 3 



= 512^ (* " Sin ^ (|Mg( ^ )|2 + |Eg|2) d ^ dcos ^ 

(3.68) 

where the definition (I2.5ip was inserted for |p+| and the relations 

= p 2 = m 2 P — 2mpE 1 = m 2 P 1 ] , (3.69) 

\ m P J 

yfs^&yfs^ = ^ds wn = — mp di? 7 (3.70) 

were used in the last line. Note that the minus sign in Eq. (I3.70P cancels against a minus sign 
resulting under a switch of the upper and lower integration limits of the dE 1 integration, 
since _E™ ax = E^(s^ n ) and vice versa, see Eq. (13.691) . Also note that equation (13.681) is 
exactly the result (6) of Ref. 
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3.5.3 The form factors for P — > 7r + 7r 7 

We have marked the form factors of Ref. [72] for the decay P — > 7r + 7r _ 7 (with a real photon) 
by a subindex G, i.e. Mq and Eq. They differ in normalization and structure from the 
corresponding form factors M and E± of Ref. [23] for the decay P —> 7i + ir~e + e~ . Let 
us first focus on the magnetic form factor Mq. Comparing the corresponding amplitudes 
fl3T62|) with fl3T80|) (i.e. Eq. (2) of Ref. [72] with Eq. (1) of Ref. [23]) we can read off the 
following relation between Mq and M: 

M G (s nn ) = m 3 P M(s nn ,k 2 = 0). (3.71) 

Here M(s n7T , k 2 = 0) contains the information of the decaying particle and the vector meson 
dominance input: 

M x VMD{s nn ) (3.72) 

if P = 7T°; 

x - 2v/2^ sin 6 mi ^j if P = 77; (3.73) 

B + 2V2^cos0 roia! ) iiP = rf. 

This means that Mg(s„) is dimensionless, whereas M(s 7r7r , 0) has the dimension [mass] -3 . 
The vector meson dominance factor contains the contributions from the terms Cpppa, 
Cpvv and also a small contribution from the Cpppy , which actually vanishes for the 
hidden gauge case. This is shown in Figures I3.15[ 13.161 and 13.171 



M(s n7r , k 2 = 0) 



with 



M 



{ 8tt^3 ^3 V h 



sin 61 





=W\A, 



Figure 3.15: direct contribu- 
tion CpppA 



Figure 3.16: direct contribu- 
tion Cpvv 



Figure 3.17: full VMD term 
£>pppv 



Therefore the vector meson dominance factor is given by 
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Note that neither form factor is normalized to one for the special cases s n7T — >• 4m 2 or 

0. The vector meson dominance factor instead is normalized to one. 
The case of the electric form factor Eq is more complicated since it differs also in structure 
(and not only in normalization) from E± from [23] . 

0. To leading order the form factor Eq could simply be put to zero. 

1. A model with a intermediate CP- violating decay P — > tt + 7t~ and a subsequent 
Bremsstrahlung corresponds to the following structure of the electric form factor 

e g my- 



7V7T 

16 e g^ 

71 



(p+ ■ k){p- ■ k) (jn\ - s^) 2 (1 - (31 cos 2 8 W ) 



^(l -««»»«.) Bf {l - (l - ^i;)-^} 
= £fe(£,,9,). (3.75) 

This model is very unlikely to lead to measurable results, since theory and phe- 
nomenology predict that the g^+w- coupling constant is very tiny: \g vn + n -\ is es- 
timated to be less than 2.6 x 10 -16 in the Standard Model (via the CKM phase), 
less than 2 x 10~ 10 under the presence of a strong 9 term in QCD, and less than 
5 x 10 -11 in spontaneous CP violating models (with more than one Higgs particle), 
see Refs. [721 [23] and references therein. 

In Ref. [72] a form factor was constructed in terms of a 4-quark-operator O. This 
operator is assumed to be unconventional, which means that it is not constrained by 
known physics. Especially, it should not contribute directly to the decay P — > 7r + 7r~ 
and the well studied K° decays. It also should be a flavor- conserving CP violating 
four-fermion operator with explicit ss quark content, such that constraints from the 
empirical bounds on the electric dipole moment of the neutron are excluded as well. 
The following short-range operator does have these features: 

O = -^-Gsiuaulhip - kYs Wfu . (3.76) 
mp 

Here G is a free, dimensionless, 'natural' model-coefficient. The latter means that 
it is of order G < 0(1). It parameterizes the strength of the operator. Geng et al. 
([72]) calculated the contribution under the assumption that the production of the 
photon results from the strangeness containing part, while the u- and (i-quark part 
is responsible for the dipion-production from the vacuum. This means: 

(P\ O |7T + 7r"7) ~ —G (P\ sia^jsip - k) u s | 7 ) (0| u^u \tt + tj-) (3.77) 

Tfl p 
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P 




Figure 3.18: Structure-dependent contribution of the unconventional operator of Ref. [72] to 

P —7- 7T + 7T~7 . 

which is visualized in Figure 13.181 

The authors of Ref [72] introduced for the P — > 7 transition the form factor F(s): 

(P\ sia^ip - kfs | 7 ) = ie[e^k ■ p) - (e • p)fc M ] (3.78) 

It can be calculated in the framework of quark models (see [201 HI] an d [ZSj)- In Ref. 
[72] it was suggested to use F(s n7T ) ~ F(0) ps 0.19. 
Assuming this value we find 

£ G (s^) ~ 2eF{s nw )G = 2eF (m 2 v - 2m v E^ G = E G (E 7 ) . (3.79) 
In our calculations we will use the form factor given in item 2. 

3.5.4 The reaction P -> tt+tt 

We will now concentrate on the decay P — > tt + 7t~PP. We start out with the invariant 
decay amplitude given in Eq. (1) of Ref. [23] : 

A(P ^ n+n-pr) = ^ e u(k- , s-)-fv(k+, s+^ Me^k^jl + E+ V % + E_p^j . (3.80) 

=jM(fc_,fc+) 

Note that the polarizations appearing in the process P —> 7r + 7r~7 are now replaced by a 
current into a lepton pair. Furthermore, we will use the following 'projection' tensor con- 
structed from the bilinear combination of the current j M (fc-, s_; k + , s + ) under a summation 
over the final spins s_ and s + of the outgoing PP pair, which can be found in Appendix 
A5: 
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If the Lorentz-indices /x, // are space-like (as it is the case for the on-shell photon, = 
(0,e*j_)), then C M/i /(/c_, k + ) gives a positive contribution. 

Note that (k + — k~)^ = for on-shell l + and l~ . Thus the spin-summed squared 
amplitude reads 



\A\i(P 7T+7T-/+/-) = jj^O»»' (Me^vXvt + E + p% + E_pt) 



2e 2 

(k 2 y 



\M\ 2 e^e^ a ,^k 2 ){-g^)ky%jlk v ' P %p 



- W e MV ^ (kl-kt) k v p%pt (k+ - kt ) k v ' p%p^ 

{ME* + + E + M*}e^(kt-k»)k»p a y_ {k+-k-,)pi 
{ME*_ + E„M*}e^(k1-k»_)k»p a y_ (k+-k-,)p^ 

\E + \*[ky + - (k ■ P+ r + ((k+ - k-) ■ p + n 

\E-\ 2 [k 2 p 2 _ - (k -p-) 2 + {(k+ -k-)- p^) 2 } 
(E* + E_ + E + E*_)[k 2 ( P+ ■ p_) - (k ■ p + )(k ■ p_) 



(3.82) 



Inserting p_ = P — p + — k into the antisymmetric products in (I3.82p and making use of 
the total antisymmetric tensor, we get 



\a\ 2 (p ^ n + n~rr 



2e 2 



(k 2 y 



\M\ 2 [e^e^ al ^k 2 ){-g^)k v plP p k v 'piP p ' 

- e^ aP e M w/j' (kl-kt) k v p%P\ki - F' ) k v 'pi PP' 

- Re{MEX}e„ v0l p{e + -kt)k v p%P^ (k+ - jti + P*) 
+ Re{MEl}e^ a p(k»-k»)k»plPP(k+-k;,)(pi -j?L - P*) 

- \E + \ 2 [k 2 p 2 + - (k ■ p + f + ((*+ - fc_) • P+ ) 2 } 

- \E_\ 2 [k 2 p 2 _ - (k ■ p^f + ((*+ - fc_) • p_) 2 ] 

- (ElE^ + E + El)[k 2 ( P+ ■ p.) - (k ■ p + )(k • p_) 

+ {k + -k_)-p + {k + -k_)-pA. (3.83) 



In the rest frame of the P meson, where P^ = {Py = m P 6^°, the spin-summed squared 
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amplitude reads 



\A\ 2 (P -»■ 7T+7T-/ + /- 



2e 2 



2^2 



ijk £ i'j'k' 



(k l + -k l _)k^ + (k l + -k l _)k 3 % 



+ Re{M(E* + )}m P e ij ' fc (fc|, -fcL)^ {K'-K')^' 



+ Re{M(E; - £;:)}m P £ yfc (^-A;l)A;^ (fcJ-^)(P 

- W + - (^-P + ) 2 + ((fef - *_) -P + ) 2 ] 

- \E_\ 2 [k 2 p 2 _ - (k ■ p„) 2 + ((k + - k.) ■ v~f] 

- (E* + E_ + E + E*_)[k 2 {p + -p„) - (k-p + )(k-pJ) 

+ (k+ - fc_) • p+(k + - ■ p_] 



p: 



(3.84) 



where we used — g M ' = 5 n> in the first term. 



Since the scalar product over four-vectors is invariant in any frame, we can drop the tilde 
in the expressions 



(k+-k-) ■ (p+ -p_) = (k+-k-) ■ (p+-p-) 



(kt-kZ)P* = (*+-*_) ■ P = {k + -k. 



P 



Formulated in terms of 3-vectors, the spin-summed squared amplitude reads then 



2e 2 



(k 2 y 



\M\ 2 m 2 P \k\ 2 (k 2 )(p + ± ) 2 - ((k+ - k][) • z A p+)' 



+ Re{M(E* + - E*_)}m P \k\ (k+-kj) • z A p+ (fc+-Jfc_) 
+ Re{M(E* + + E*_)}m 2 P \k\ (k±-kjj -zAp+(fc + - k. 

- \E + \ 2 [k 2 p 2 + - (k- p + ) 2 + ((k + - • p + f] 

- \E„\ 2 [k 2 p 2 _ - (k ■ p_) 2 + ((k + - fc_) • p_) 2 ] 

- (E* + E„ + E + El)[k 2 (p + ■ p_) - (k ■ p + ){k ■ p_) 

+ {k + -kJ)- P+ {k + -kJ)-pj\. 



(P+-P-) 
■Pv 



After (i) the relations flA.8[) and flA.9j) are applied, (ii) the scalar products (k + — Af~) • (p + — 
p~) and (k + — k~)-P are simplified with the help of (IA.25|) and (1A.26[) . respectively, and (iii) 
the relations (IA.30p . 0A.3ip and ( 1A.34j) are inserted in the electric terms, the spin-summed 
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squared amplitude reads 



|^| 2 (P -»• 7T + 7T-/ + / 

2e 2 



|M(s^)| 2 m||k| 2 |p* | 2 sin 2 ^ 



(k 2 y 



(k ) - 4|ki| sin ^sin> 



+ Re{M(£;* -^I)}m P /3 7r /9 z |k||p* ||kl| sin ^ sin 0, 
|-(mp-s„-S|() cos^,!- cos sin <^ + 2 v /s 7r7rV /s» sin 9 n sin 6 1 ; sin <p cos <^ 

+ Re{M(E; + E*_ ) }m P \k\ \ p* + \ | kl | ft A 1/2 (s„ , m 2 ,, s™) sin W sin 0j cos 9 X sin </? 



1*4 



16 



A 1 / 2 (s 7r7r , ffjp, fc 2 ) — (m 2 p — s n7T — /c 2 ) (3 n cos 



(l-/3?cos 2 0,) 



+ sin 2 6^ (l - /3 2 sin 2 9 L cos 2 <p) 

A 1 ^ 2 (s 7r7r , mp, k 2 ) — {m 2 P — — k 2 )(3 n cos 9^ 
1 



's^finy/sufii sin sin cos 9\ cos <£> 



16 



A 1//2 (s 7r7r , mp, /c 2 ) + (m 2 p — s nn — k 2 ) (3 n cos 9 n 



(l-PfcoffOt) 



+ \s nn su/3* sin 2 9 n (l- (3 2 sin 2 9 X cos 2 </?) 



+ 1 



A 1 ^ 2 (s vr7r , m 2 P , k 2 ) + (mp — s nn — k 2 )f3 n cos 9 V 



/ S n n(3n\/Suf3'l Sill 9 n Sm 9\ COS 9\ COS 



(E* + E_ + ) ^A(s^, m 2 , k 2 ) (l - dl cos 2 6 A (l - ft 2 cos 2 9 t ) 

— jS^suP 2 I 1 — Pf sin 2 6*^ sin 2 9i cos 2 <^ — /5f cos 2 9^ cos 2 ^ 



— zV ' s ™y/~siiftiPl ! (% — s^tt — A; 2 ) sin 6^ cos sin 6^ cos 6^ cos </> 



(3.85) 



We also used 

cos <£ p+ (- sin ) - sin (- cos y? fc _ 



sm 



(V3 P+ - ip k _) = + sirup 



and 



Mp+AkD = (pt) x (*i) w - (pt) w (*ir 

= Ip+lsin^Tf |k*|sin#/ (cos sin 99^- — sin <p w + cos 



IpII 



|kH 



— I p + 1 sin 9-n I k* I sin 9\ sin 



(3.86) 



where tp is the angle between p+ = — p ± and k+ = — k ± in any frame (especially also in 
the 77 rest frame). Since pjj", pjT, kjj" and kjj~ are parallel or antiparallel to each other and 
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to the positive ^-direction, tp is also the angle between the l + l plane and the ir + n plane 
(see the text below Eq. (1) of Ref.[23]). 

Expressed in terms of the standard variables the squared matrix element for the decay 
P — > 7T + 7r~Z + / _ reads: 



s u ,9 w ,9 h p) 

e 2 \ 

= 8(fc 2 ) 2 | l^( s ™> s ")l 2 ^( m P> s '"" s ") [1 - 0? sin 2 9 t sin 2 p] s wn /3 2 sin 2 9 n 

+ 4Re{M(s„, s n ){E* + - E*_)} \ 1/2 (m% s^, s u )Pf ftV^V^i 

x f — ~ {m 2 P — s n7T — su) sin 9^ cos 9 W cos 6i sin 9i sin + y/s nn sn sin 2 0^ sin 2 0; sin p> cos 99) 
+ 2Re{M(s 7r7r , s«) (£7+ + ) } ^/s^s u f3 2 (3 n \(m 2 P , s nn , s u ) sm9 n sin^ cos9 t sirup 

+ 4s 7r7r s // /3 2 sin 2 9 n (l - /3 2 sin 2 6 1 , cos 2 y?) 
=F 4 [A 1/2 (s 7r7r , m F ] y/s^(3 wy /sTi(3 2 sin 6^ sin 6 1 / cos Z cos p 

+ 2Re[£;£_] A(s 7r7r , m 2 P , k 2 ) (l - /3 2 cos 2 W ) (l - ^ cos 2 0,) 

— As^suP 2 (l — f3f sin 2 0„. sin 2 9i cos 2 <£> — /3 2 cos 2 0^ cos 2 0;) 

— 4y/s W7T yfsu B 2 B 2 (mp — s wn — k 2 ) sin 9 n cos 0„. sin 0; cos 9\ cos <^ 



(3.87) 



3.5.5 Decay rate 



We will now insert the expression of the squared matrix element ( I3.87P into the formula of 
the decay rate (1A.46|) . 



Magnetic term 

We will start with the magnetic term. As we will see in the next Chapter, it is the leading 
contribution. The final expression (with k 2 = su) is: 

dr(|M| 2 ) = j^- 2 \M\ 2 2 * 12 6 m 2 P K\ \pl\ 3 \k\ 3 k 2 [l-(3? sin 2 0, sin 2 <p] sin% ( Smr )"3 

x ds n7T dcos 9i dcos 9^ dp dk 2 
\M\ 2 I J 3i (l - \f3 2 ) s 7r7r /3^A5(mp,s 7r7r ,A; 2 )sin 2 6' 7r ds 7r7r dcos6' 7r dA; 2 



fc 2 ' 1 m P 2 16 vr 5 11 \ 3 
e 2 1 3 3 - B 2 

— \M\ 2 3 3 g g^A^mp, s™, — ds^r ds H . (3.88) 
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This expression exactly agrees with Eq. (4) of Ref. |21j . 



Mixed term 



With the same definitions as before we can calculate the decay rate of the mixed term, 
weighted by sign(sin<y?cos<£>). Note that only the second part of the mixed term (oc 
~Re[M(E+ — -E*)]) will give a contribution. The other two terms will be equal to zero 
for the following reasons: (i) they vanish if we integrate over the angle 6 n (and also over 
9i in case of the last term), because of the cos 9^ sin 8 n dependence; (ii) the main reason 
- valid also for the case of generalized form factors - is the integration over the angle ip: 
with the weight sign(sin<y9cos</9) this integration reduces all of the above mentioned terms 
to zero, unless there is the structure sinyjcost^. 

Note that this is actually the CP-violating term because the remaining part is proportional 
to Re[M(E+ — EZ.)], which is in fact CP- violating. The latter term would vanish, as in the 
P — y 7r + 7r - 7-decay, if the norm | sinyj cos ip\ were not taken. This term can be measured by 
analyzing the forward-backward asymmetry of the angle </>. If is not zero in average, there 
does exist a term proportional to Re[M(E+ — E*_ )] and so CP- violation can be measured. 



The decay rate for the mixed term is given as 



dr(Re((£ + - £_)Af*)sign(sin<pcosy>)) 

:Re[M(£;-^)]| P :;| 3 |k| 2 |k _| 3 



si m2 3V /%2i%e 
x sin 2 ^ sin 2 ^ | sin ip cos (p\ (s^) -5 ds W7T dcos 9i dcos 9 n dip dsn 

\(mp, s^r, su)s^l sin 2 9 n 0fds nn dcos ff ds u 



e 2 Re[M(E* + -El)]^ 



3 • 2 13 vr 6 mp s u 
e 2 Re[M(E* + -El)}X(m 2 P ,s 7T7T ,s ll ) 
3 ■ 2 13 7T 6 mp s u 



s n7T f3l sin 2 ^ ds nn dcos 9^ dsu . (3.89) 



In the last step the approximation /?/ ~ 1 was inserted. This can be justified for the case 
that the leptons are electrons. 

The CP-violating observable Aqp which can be measured by the experiments is given by 
the mixed term normalized to the total decay width (which is given to leading order by 
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the integral of Eq. 02 



.4 



Jo 



2k 



dT(rj — > tx + ti e + e~ 



-d0sign(sin0cosi 



CP 



2-rr dT(r] — > n + ix e + e~ 



d<f) 



3 ■ 2 13 7r 6 mp T(?7 7r+7T-e+e 
A(mp, s^, s« 



X 



su 



s,™/3* sin 2 ^ /3 Z Re [M ( ST — E*_ ) ] ds^ dcos ^ ds u . 



(3.90) 



Note that our result for Aq? agrees with Eq. (3) of Ref. [23] when the mentioned approxi- 
mation 0i f=s 1 is inserted. 



Electric terms 

The decay rates of the electric terms can be calculated in the same way as the ones of the 
magnetic and mixed terms before. 

In the following, we will separately construct the decay rates of the squared electric terms 
which are proportional to |-E±| 2 and the decay rate of the mixed electric terms proportional 
to (E* + E_ + E + E1) = 2Re[E* + E_}: 



dT(\E ± \ 2 ) 
e 2 \E + \ 2 



A 1/2 i 



k 2 ) =F (m 2 P - s n7r - s u ) (3 W cos 6 V 



1 -/3 2 cos 2 z 



+ As n7T su(3l sin 2 6 lx [l — f3 2 sin 2 6i cos 2 ip ) + 0(sin 6 e cos Q\ cos ip) 
1 



x 



mp ■ 2 15 ■ 7r 6 



/3 7 rAA 1/2 (s 7r7r ,mp, s u ) ds n7T dcos 6», dcos 6 n dipds a 



e 2 \E±\ 2 



3-/3/ 



s 2 ■ mp ■ 3 • 2 15 • 7T 5 

: [A 1/2 (s„,mp 



s/0 =F (mp-s^-sji) cos 9 n ] 2 + As^suftl sin2 
x /3 7r /3,A 1/2 (s 7r7r , mp, s M ) ds^ dcos 6^ ds H . (3.91) 



The mixed term has a slightly different form compared to the previous expression and 



3.5. P -> 7T+7T-7 AND P tt+tt^L+L- 



49 



reads: 



A(s 7r7r , m P , s H ) 1 - cos 6^ • (1 - (9, cos 9 t ) 



dT((ElE^ + E + E*_)) 
e 2 {E\E_ + E + E*_ 

— 4s 7r7r s«/3 2 ^1 — /3f cos 2 6*^ cos 2 9 { — f3f sin 2 0^ sin 2 9i cos 2 </? ) + 0(sin 9i cos cos <p) 
1 



x 



mp ■ 2 15 • 7r 6 



ArA-^^TrTr, m p , S») ds ff7r dcOS dcOS 6^ d<£ dsjj 



e 2 Re[E;E_] 3 - # 



A(v,mp, s«) I l-^cos 2 ^ ) - 4s n7r siif3„ 



(k 2 ) 2 ■ mp ■ 3 • 2 14 • vr 5 2 

x (3 n (3i\ 1/2 (s nir , m|, s M ) ds^ dcos 9 n ds a . 

3.5.6 The form factors for the decay 77 — > tt + tt 

Here we discuss the magnetic and electric form factors M and E for the decay P — > 
7r + 7r _ / + / _ as defined in Ref. [23] and as used in the last Section. The form factors are 
modeled according to the one of the P — > 7r + 7r~7-decay, but as mentioned before, they 
differ in the powers of masses. 

We start again with the magnetic form factor M(s 7rvr , su). This is again: 

M(s^, s u ) = Mx VMD(s 7rw , s„) (3.92) 
where VMD is the pertinent vector meson dominance factor and Ai is given by: 

6 if P = 7T°- 

8vr 2 /3 ' 
-^^-y=^cos9 mix -2V2^sin9 mix ^ if P = 77; (3.93) 

g^3^= sin9 mix + 2V2j^ cos 9 mix ^j if P = 7/ 



M 



with the pion decay constant f n ~ 92.4 MeV, the octet pseudoscalar decay constant /§ ~ 
1.3 f n , the singlet pseudoscalar decay constant / ~ 1-04/^ and the 77-77' mixing angle 
#mix ~ -20°, see Ref. [32]. 

The VMD form factor again contains contributions from the terms Cpppa, £pvv and 
Cpppv as it was already shown in figure 13.151 figure 13.161 and figure 13.171 Especially 
VMD(s n7T , Su) differs from VMD(s nn ) of Ref. [23J by an additional form factor related to 
this part of the off-shell photon (decaying into the I + 1— pair) that according to Ref. 
does not directly couple to the Pti + it^ complex, see Ref. [21]: 
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Therefore the vector meson dominance factor has the following form: 

3 , 3 . s mi 



VMD^s^su) = 1 - -{ci - c 2 + c 3 ) + -{ci - c 2 - c 3 ) 2 



rriy m 2 



+ -C3^ — —5 — ^7 r. (3.95) 

2 - s a - im v T(su) my - s wn - im v T(s wn ) 



By adjusting the values of the Cj-parameters, we can switch between the various VMD 
models. Note that for the hidden gauge case the Cpppy term vanishes. 
The electric form factors E± which contribute via the combination R,e[M(E+ —E*L)\ to the 
squared CP-breaking amplitude are model-dependent. 

0) To leading order, E± can be put to zero. 
In Refs. [72| 123] two models for E± can be found: 

1) The first model consists of induced Bremsstrahlung of an 7r + 7r~ intermediate state 
which violates CP symmetry (see Eqs.(7)-(14) of Ref. [23]): 



k 2 + 2k ■ p± 

4eg P7r+7r - 



(3.96) 



m 2 p + q 2 -s n7T =F /3 7r A 1 /2( m 2 j) Snnj g 2) C0S Q n ' 

As argued in the previous Section, this model will not play a role in our discussions. 

2) The second model of Refs. [721 123] describes the CP violating P — > n + 7T~l + l~ decay 
to a short-range E\ operator, which was constructed for the P — > 7r + 7r~7-decay. 
According to Eqs. (15-17) of Ref. [23] the electric form factors E± have the following 
parameterization: 

£±(WaA) = ± eF(S 7 )G (k 2 + 2k-p T ) 

Tfl p 

& F(S-mr) G / / o \ n i /o / 9 \ \ 

(± (m P + s« - s^J + p^A 7 (m P , s^, s u ) cos ff J 



2mp 



(3.97) 



with the form factor F(s nn ) ~ -F(O) ~ 0.19 and G a free model coefficient of order 
G < 0(1). 

Note that in this case 

eF(s JTn )G 2 

Ej r (s 7T7T1 Su, 6 W ) £/_ (St,-,,-, S^, 6 w ) = 3 {fTlp + Sll — s n7T ) 

= E+{sm, s u ) - E_{s^, s u ) (3.98) 
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which is in fact independent of W , while 



|2 



= 1 6 2 i ( m P + S «- S 7T7rj + P n \{m p , S nn , Sii) COS 0* 

=F 2s^r (m 2 p+su-s nn ) (3 n \ 1/2 (m 2 P , s n7T , s u ) cos9 n 

= \E(s W7 r, su, cos 2 9 n )\ 2 ± AE 2 (s n7r ,s u )cos9 n (3.99) 

contains even and odd powers in cos^. 

In general, the difference between E + (s n7T , su, 9 n ) and E_(s n7T , su, 9 n ) is a function of only 
even powers of cos 6^, i.e. 

E + (s n7r ,su,9 n ) - E^s^su.O^) = A£ ± (s„, s u , cos 2 ff ) , (3.100) 

whereas the squared moduli of i? + and E_ have the following dependencies: 

I E + {s^, si U Bk) 1 2 = I-B^jSh, cos 2 0^)1 2 - A J E 2 (s 7r7r ,s Zi ,cos 2 7r )cos0 7r , 
|^_(s^, s«, 0^)| 2 = ^(sttt,-, cos 2 9 n )\ 2 + AE 2 (s nn , su, cos 2 9 n ) cos9 n . 

(3.101) 
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Chapter 4 



Results 



In this Chapter we will state the results for the branching ratios of pseudoscalar mesons 
P = 7T°, rj, rf for all discussed modes. The calculations have to be done numerically. In 
practice, we applied the Gaussian routine of the CERNLIB library [76] to perform the 
integration over the squared invariant masses of the outgoing particles. The decays into 
three particles, namely P — > / + /~7 and P — > 7r + 7r~7, can be handled easily, because one 
has to integrate over just one variable. In comparison, the decays into four particles, 
P — > l + l~l + l~ and P — > 7r + 7r~/ + /~, are more complicated, since one has to integrate over 
two variables which depend on each other. 

We compare the results for the different vector meson dominance models (the hidden 
gauge and the new modified one) with the experimental data and other already published 
theoretical calculations. 

The listed errors can be traced back to the different vector meson masses which occur in the 
VMD factor. As pointed out in the second Chapter, see Eq. (12.361) . the coupling constant 
g is related to the vector meson mass. So according to the different fits given in [30] various 
vector meson masses are used in our calculations to take this into account. For the hidden 
gauge model the vector meson mass varies from my = 760 MeV to my = 782 MeV. For 
the modified model we use a vector meson masses from my = 760 MeV to my = 791 MeV, 
because the relevant fits generate very high coupling constants. 

As for the modified model, there exist also a small contribution to the total uncertainty by 
the ambiguity of the two fits of Ref [5U] . This contribution will be in range of a couple of 
per cent of the total errors, so that we will not notice them in the most cases. Because we 
normalize the decays P —> / + /~7 and P — > l + l~l + l~ to the two-photon decay, there appear 
no errors from the mixing angle of r/o and r]g and the coupling constants fo and fg. They are 
indeed present when the absolute decay rates for the P —> 7r + 7r~7- and 7r + 7r~/ + /~-decays 
are given. 
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4.1 P -> Z+Z"7 



We start with the discussion of the decay of a pseudoscalar meson into a lepton-antilepton 
pair and a photon. We show the dependence of the branching ratio expressed relative to 
the decay into two photons on the invariant masses of the leptons and give the branching 
ratios expressed relative to the decay into two photons: 



BR rel (P Z + Z"7 W ) 



BR(P -> Z+Z~ 7 
BR(P 77) 



(*)^ 



(4.1) 



7r u — > e + e 7 

The dependence for the relative branching ratio of 7r° — > e + e~7 on the invariant mass of 
the electrons, 



is shown in Figure 14.11 




[MgV] 



Figure 4.1: Dependency of the differential relative branching ratio (normalized to BR(7r° — > 77), 
see Eq. (|4.ip ) of the decay 7r° — > e + e~7 on the change of the invariant mass of the electron pairs. 



The strong peak for small masses is typical for branching ratios plotted against the invariant 
mass of electrons, y/s^, as we will see in following sections. The errors are smaller than 
the width of the lind_|. Also the different models, the one without VMD, the hidden gauge, 
and the modified VMD model all give the same graph. 

In Table 14.11 the results for the branching ratios of the decay 7r° — > e + e~7 are listed for 
the cases (i) without a VMD factor, (ii) with the hidden gauge model, and (iii) with the 
modified hidden gauge model. The errors are smaller than the number of decimal places, 
so we do not show them. 



1 Here and in the following calculations based on two on/off-shell photons the differential branching 
ratios are calculated as ratios to the total 2-photon decay rate, such that the coefficients of the /o, fs, and 
the mixing angle 9 m i X drop out. The plotted values are therefore dimcnsionlcss. 
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without VMD 


hidden gauge 


modified VMD 


BR rcl (7T° -> e+e- 7 )(10- 2 ) 


1.185 


1.188 


1.187 



Table 4.1: Relative branching ratios (normalized to BR(7r° — > 77), see Eq. (|4,ip ) of the decay 
7T° — > e + e~ 7 calculated without VMD, with the hidden gauge model, and the modified VMD 
model. 

There is a small difference between the values with and without a VMD factor, but the 
difference between the two VMD models is hardly noticeable. 

In Table H~2l we list the corresponding values calculated by other groups and compare them 
to the experimental data. 





m 


m 


i 


exp. data [34] 


BR rcl (7r°^e + e-7)(10- 2 ) 


1.18 


1.18 


1.18 


1.188 ±0.034 



Table 4.2: Published theoretical values and experimental data for the relative branching ratios 
of the decay 7r° — > e + e~ 7 . 

A comparison of Table 14.11 and 14.21 confirms that our results agree with the published ones 
and are consistent with the current data [34] . 

77 — > / + /~7 

The next results which we present are the ones of the decays rj — > e + e~7 and 77 — > /i + /i~7. 
The dependences on the invariant masses of the dileptons are shown in Figure 14.21 

The curve plotted against ^fsZ e is almost the same as in the analogous 7r°-decay, except for 
the amplitude. The graph plotted against y/s^ is wider and also has a smaller amplitude. 
The errors are again smaller than the width of the lines, but the contribution of the VMD 
factor is visible now. 

The branching ratios for these decays are listed in Table 14.31 

Qualitatively, the results of the decay 77 — > e + e~7 are similar to the ones in the 7r° sector. 
The differences between the values calculated without VMD and the with VMD are still 
small, while the two VMD models generate almost the same results. 

For the decay 77 — > /i + yU^7, however, there is a clear discrepancy between the predictions 
with and without VMD. The two VMD models are still in the same range, although the 
difference has increased to 0.3% now. 

Other theoretical values and experimental data are listed in Table 14.41 
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Figure 4.2: Dependence of the differential relative branching ratio (normalized to BR(t] — > 77), 
see Eq. (|4.ip ) of the decays rj — > e + e~j and rj — > /i + /i~7 on the change of the invariant mass of 
the electron-pairs and the muon-pairs, respectively. 





without VMD 


hidden gauge 


modified VMD 


BR rcl (r? -> e + e- 7 )(10^ 2 ) 


1.619 


1.666 ±0.002 


1.662 ±0.002 


BR rcl (r/^/x+/i"7)(10- 4 ) 


5.51 


7.75 ± 0.09 


7.54 ±0.11 



Table 4.3: Relative branching ratios (normalized to BR(?] — > 77), see Eq. ()4.ip ) of the decays 
7) — > e + e~7 and r/ — > ri + fi~-y calculated without VMD, with the hidden gauge model, and the 
modified VMD model. 





m 


M 


exp. data [54] 


BR rcl (r/^ e + e- 7 )(10- 2 ) 


1.62 


1.77 


1.78 ±0.19 


BR rcl (r/^/i+/i"7)(10- 4 ) 


5.54 


7.48 


7.9 ± 1.1 



Table 4.4: Theoretical values and experimental data for the relative branching ratios of the 
decays 77 — > e + e~j and rj — > fi + fj,~^f. 



The values given in Ref. [7] agree with our calculation without a VMD factor. These 
values fit the data for the decay 77 — > e + e~7, but have no overlap with the experimental 
data of the decay 77 — > /i + /i~7. For the case of the 77 — > e + e~7-decay, the values of [9] are 
larger than ours, while both calculations are consistent for the 77 — >■ /i + /i~7-decay. Note 
that our VMD values represent the data very well, compare Table l4~3l and I4~4l 
The authors of [77] calculated also the corresponding branching ratio in the framework of 
VMD, i.e. with the hidden gauge model and without a VMD factor, respectively. The 
corresponding values are totally analogous to ours. 
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Finally, we discuss the results for the 7/ decays. The dependence of the differential branch- 
ing ratio of the decay rf — > e + e~7 looks basically the same as in the rj sector. As shown 
in Figure 14.31 there is again a very strong peak for small energies. But there is also a 
contribution from the width, a bend in the region of the vector meson mass. To visualize 
that we plotted also the relevant sector (see Figure l4~3j) . 




Figure 4.3: Dependency of the differential relative branching ratio (normalized to BR(r/' — > 77), 
see Eq. (|4.ip ) of the decay rf — > e + e~7 on the change of the invariant mass of the electron-pair 
in different energy regions. 



In the decay rf — > 7 the contribution of the width is very large. It is shown also in 
Figure l4~4l 
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Figure 4.4: Dependency of the differential relative branching ratio (normalized to BR(r/' — > 77), 
see Eq. (|4.1|) ) of the decay rf —> ^ + /i _ 7 on the change of the invariant mass of the muon pair. 
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The branching ratios are listed in Table 14.51 





without VMD 


hidden gauge 


modified VMD 


exp. data [34] 


BR rcl (r/ ^e+e- 7 )(10- 2 ) 


1.79 


2.10 ±0.02 


2.06 ±0.02 


< 4.5 


BR rcl (r/ ^/i+u"7)(i(T 3 ) 


1.72 


4.45 ±0.15 


4.11 ±0.18 


4.9 ± 1.2 



Table 4.5: Relative branching ratios (normalized to BR(?/ — > 77), see Eq. (|4.ip ) of the decays 
7/ — > e + e~7 and rj — > /i + /U _ 7 calculated without VMD, with the hidden gauge model, the 
modified VMD model, and the experimental data. 



The trend of the r\ decay channels is repeated here. The values of the various VMD models 
hardly differ, while the values calculated without VMD are clearly smaller. 
All results are compatible with the upper experimental limit of the 7/ — > e + e~7 decay. For 
the case of the decay rj' — > /i + yU~7, the value calculated without a VMD factor is out of 
the experimental range. Both VMD results are smaller than the experimental data, and 
only the hidden gauge model has a small overlap. This may be significant, but notice that 
the experimental accuracy in the 7/ sector is also not very high. 



Summary 

In general, the theoretical values represent the data very well. The contribution of the 
VMD factor is very small for the decays P — > e + e~7 and increases with the mass of 
the decaying particle. One can see that for the decays 77/77' — > fi + fi~ , -f a vector meson 
dominance model is needed for an accurate description of the data, but no preference of 
any model can be given. 



4.2 P -> /+/"/+/- 

We will now present the results for the rates of the decays of a pseudoscalar meson, 
P = n°,rj,ri', into four leptons. 

The direct term and the crossed term give the same result under integration, so we just 
need to calculate one of them. The interference term is more complicated to calculate, 
because there appear functions of the three respective angles, namely 6*12, 6*34 and <fi see 
Eq. f)3.44p . fIB.lj) . flB.2j) and flB.5j) . in the denominator. Thus the angle integrations are 



no longer trivial. So instead of a two-dimensional integration the task now is to perform a 
Gaussian integration over five variables. 
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7T 



e T e e^e 



We start with the four lepton decay of the 7r°. Here only the decay into e + e~e + e~ is 
possible. The dependence of the branching ratio on the invariant mass of the electron (y/s^) 
is shown in Figure 14.51 Note that the x-axis is of logarithmic scale and the interference 
term is scaled up by a factor of —100. 
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Figure 4.5: Dependency of the relative branching ratio (normalized to BR(-7r° — > 77), see Eq. 
(|4.ip ) of the decay tt° — > e + e~e + e~ on the change of the invariant mass of the electrons. 



For the direct (crossed) term we can see a peak at very small ^fs^. The interference term 
in contrast is negative and wider. The different VMD models give almost the same curves 
for the direct (crossed) and the interference term. So we just plotted one curve, respec- 
tively. The errors are smaller than the width of the lines. 

The calculated values of the branching ratios are given in Table I4.6I The errors are again 
very small so that we do not show them. 



7T° e+e-e + e~(10- 5 ) 


without VMD 


hidden gauge 


modified VMD 




3.456 


3.469 


3.468 




-0.036 


-0.037 


-0.036 


nprel 
DIX total 


3.420 


3.432 


3.431 



Table 4.6: Relative branching ratios (normalized to BR(7r° — > 77), see Eq. (|4.ip ) of the decay 
7T° — > e + e~ e + e~ calculated without VMD, with the hidden gauge model, and with the modified 
VMD model for the direct and crossed term BR^, the interference term BR^ 1 , and the total 
one BRlf tal = BR\% + BR#. 



As one can see the contribution of the interference term (BR^ ) is of the order of one per 
cent compared with the leading direct and crossed term (BR^ 2 ). Note that the inter- 
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ference term has a negative sign. The difference between the VMD models is very small 
for this decay. We display four decimal places just to see a difference in the values of 
the interference term. Moreover the errors are small. This is compatible with the small 
contribution of the VMD. 

We can compare our results with other theoretical calculations and the experimental data, 
see Table EZJ 



7T° ->■ e+e-e+e-(10- 5 ) 


m 


m 


i 


i 


Data [31] 


bK l+2 




3.46 


3.456 






£>K, 12 




-0.18 


-0.036 






nnrel 
D5x total 


3.47 


3.28 


3.42 


3.29 


3.38 ±0.16 



Table 4.7: Theoretical calculations and experimental data for the relative branching ratios of 
the decay tt° — > e + e~e + e~. 

For the case of the total rate we agree with all values of the recent references, namely [6], 
[7j,[B] and pj. Note that the comparison has to be done with the first column of Table 
14.61 because none of them used a vector meson dominance model. With respect to the 
interference term, we totally agree with the values given in j8]. In agreement with this 
reference we also differ by a factor of 5 compared with the result of Ref. [7]. All of the 
predictions for the total branching ratio are consistent with the experimental data. 



rj^l+l-l+l- 

In the case of the rj, all three decays into two e + e~-pairs, two /i + /i _ -pairs and an e + e~- 
plus a /i + /i~-pair are possible. The dependence of the differential branching ratio on the 
invariant masses can be seen in Figures 14.71 and 14.61 

The decay into two e + e~-pairs is similar to the analogous decay of the tt°. The errors are 
again smaller than the width of the line, but we can see the difference between the used 
models. The curve without VMD is smaller than the ones calculated with VMD models. 
For the direct and crossed term we can also see a small difference between the hidden gauge 
and the modified VMD model. 

The graph for the decay into two /i + /i~-pairs is very different. We do not need a logarithmic 
scale to present the values appropriately because the direct and crossed terms are much 
wider. Also the amplitude is smaller. The difference between the used models is obvious, 
even for the interference term. For the first time we can see the errors in the direct and 
crossed terms. 

The differential ratio for the decay t] — > fi + fi~e + e~ plotted versus ^fs^ e has basically the 
same structure than the decay into two e + e~-pairs. Here only the amplitude is smaller. 
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Figure 4.6: Dependency of the differential relative branching ratio (normalized to BR(r] — > 77), 
see Eq. (|4.ip ) of the decays ry — > e + e~e + e~ and r/ — > on the change of the invariant 

mass of the electrons and the muons, respectively. In the left graph a) is the curve for the direct 
and crossed term calculated without a VMD factor, b) the one calculated with a VMD factor, c) 
is the curve for the interference term calculated without VMD factor, and d) the respective one 
calculated with a VMD factor. The values for the interference terms were scaled up by a factor 
—1000. The fluctuations in this curves are due to the imprecision of the Gaussian integrations. 
In the right graph a), b) and c) are the curves for the direct and crossed terms calculated without 
a VMD factor, the hidden gauge model, and the modified model, respectively, d) and e) are the 
curves of the interference term calculated without a VMD factor and with a VMD factor. The 
curves of the interference term were scaled up by a factor —5. 
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Figure 4.7: Dependency of the differential relative branching ratio (normalized to BR(r/ — > 77), 
see Eq. (|4.1|) ) of the decay rj — > fi + fj,~e + e~ on the change of the invariant mass of the electrons 
and the muons, respectively. 



On the other hand, the differential ratio plotted versus -Js^ is similar to the decay into 
two yU + /i~-pairs, with bigger and wider amplitude. 

Now we present our results of the branching ratios for these three decays, calculated with- 
out VMD, with the hidden gauge, and with the modified model, and compare them with 
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other theoretical values as well as the experimental data (Table UTS]) . 





without VMD 


hidden gauge 


modified VMD 


m 


i 


Data [31] 


7] e + e"e + e _ (10~ 5 ) 


BK 1+2 


6.497 


6.848 ±0.012 


6.817 ±0.015 


6.50 


6.26 


< 17.5 




-0.034 


-0.048 ±0.001 


-0.047 ±0.001 


-0.36 


Dnrel 
Dsx total 


6.463 


6.800 ±0.013 


6.770 ±0.016 


6.14 


v -> /^ + / 




BR rel 


3.996 


5.616 ±0.063 


5.465 ± 0.079 


1.99 


1.48 


< 40.5 


?7 — >■ (10 9 ) 


BK 1+2 


6.560 


10.031 ±0.129 


9.708 ±0.162 


6.73 


4.27 


< 9.1 • 10~ 4 


™ 12 


-0.049 


-0.078 ±0.001 


-0.074 ±0.001 


-0.50 


Dprel 


6.511 


9.953 ±0.053 


9.634 ±0.163 


6.23 



Table 4.8: Relative branching ratios (normalized to BR(?7 — > 77), see Eq. ()4.ip ) of the decay 
77 — > l + l~l + l~ with different VMD models, other theoretical values [7J, [9] and experimental data. 

M 



The relative branching ratio of the 77 — > /! + /i~/! + /i - -decay is much smaller than the branch- 
ing ratios of the other 4-lepton decays. This is caused by the very small phase space in 
this reaction. The interference term is again of the order of one per cent. 
In Ref. [77] the decays 77 — > e + e~e + e~ and 77 — > /i + /i~e + e~ were calculated with VMD 
models and without. These values totally agree with ours. 

Our values of the decays r\ — > e + e~e + e~ and r] — > fi + fi~fi + fi~ calculated without a VMD 
form factor approximately agree with the results given in [7J. The difference may be caused 
by the improved input data that we used. We also agree with [9] for the rj — > e + e~e + e~- 
decay while their values for 77 — > fi + fi~fi + fi~ are much smaller. We disagree again for the 
values of the interference term. Our interference term is here about 10 times smaller than 
the values given in [7J. Ref. [9] did not calculate an interference term. We also disagree 
in the total branching ratio of the decay 77 — > fi + fi~e + e~, but the factor 2 of [7] may be 
due to a typographical error in their calculations. This was already pointed out by the 
authors of [8], when they compared the respective results of the decay K L — > /i + /i~e + e~. 
It is clearly visible that Ref. [5] gives even smaller values than [7J. 

The vector meson dominance factor now gives a big contribution. Especially for the 
r] — > fi + fi~e + e~- and 77 — > yU + /i~yU + /i~-decays this factor is essential. While the results 
of the two VMD models are still very close to each other, the one without VMD is dis- 
tinctly smaller. Unfortunately the existing experimental data only give upper bounds 
which all of the models can meet. More precise measurements are needed to falsify the 
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VMD model. 



Also for the rj' case all three decay modes (e + e~e + e~, and /i + /i~e + e~) are 

possible. The dependence of the branching ratios on the invariant mass is almost the 
same as in the analogous ?7-decays. In general, the amplitudes in the ^'-decays are greater. 
Also the differences between the VMD models are bigger. This holds also for the errors. 
Differences can be seen in the curve of the decay rj' — > fi + fi~e + e~ , where the width gives 
a very strong contribution to the branching ratio in the region of the vector meson mass 
when plotted against the invariant mass of the muons This can be seen in Figure 




200 300 400 500 600 700 800 



Figure 4.8: Dependency of the differential relative branching ratio (normalized to BR(r/' — > 77), 
see Eq. (|4.1[) ) of the decay rj — > / Li + /x _ e + e~ on the change of the invariant mass of the muon 
pair. 

Plotted against the invariant mass of the electrons one can not see the contribution of the 
width. Because of the high mass of the muons the energy of the electron pair does not get 
in the region of the vector meson mass. 

In the curve of the decay rj — > e + e~e + e~ the contribution of the width gives a small bend 
in the region of the vector meson mass but this is hardly noticeable because the values in 
this area are very small. Therefore we plotted the respective range in Figure 14.91 
For the decay rj — > ji + ji~ ji + ji~ again the width does not play a role, because the energy of 
the muons is lower than the vector meson mass. Therefore the dependence of the differential 
branching ratio on the invariant mass is similar to the one in the decay 77 — > fi + fi~ fi + fi~ . 
The branching ratios corresponding to the 7/ decays are given in Table 14.91 

The contribution of the interference term is again of the same range. Because of the larger 
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Figure 4.9: Dependency of the differential relative branching ratio (normalized to BR(?/ — > 77), 



see Eq. (|4.1|) ) of the decay rj 



e 1 e e 1 e 



on the change of the invariant mass of the electrons. 





without VMD 


hidden gauge 


modified VMD 


77' e + e-e + e-(10" 4 ) 


BK 1+2 


0.7972 


1.0445 ±0.0136 


1.0148 ±0.0170 


BK 1+2 


-0.0033 


-0.0110 ±0.0004 


-0.0104 ±0.0003 


pnrcl 


0.7939 


1.0335 ±0.0140 


1.0044 ±0.0173 


rf ->■ / u + /i"e + e~(10~ 5 ) 


BR rel 


1.462 


3.739 ±0.132 


3.458 ±0.160 


77' /x+/x"^ + /^"(10- 6 ) 




0.443 


1.205 ±0.047 


1.119 ±0.056 


BK 1+2 


-0.050 


-0.182 ±0.003 


-0.172 ±0.011 


pnrcl 


0.393 


1.023 ±0.050 


0.947 ±0.067 



Table 4.9: Relative branching ratios (normalized to BR(?/ 
rf -> l + l~l + l~ with different VMD models. 



77), see Eq. (|4.ip ) of the decay 



rf mass and thereby a bigger phase space, the difference between the rf — > fi + fi~ fi + fi~ 
decay mode and the other modes is smaller than for the 77-decays. We have again a huge 
effect of the vector meson dominance factor for the rf — > e + e~ fi + fi~ case, and even more so 
for the rf — > fi + fi~fi + fi~ -decay. However, the theoretical and experimental status is much 
worse than for the 7r° and rj cases, since no data exist. 
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Summary 

In summary, the effect of VMD again increases with the value of the mass of the decaying 
particle and the lepton mass. It is very small for the case of the decay into four electrons 
and increases roughly by a factor of 1.5 or 2.5, respectively, for the cases of the decay into 
four muons. The effect of the modified VMD model is noticeable, but not very large. A 
test whether a VMD model is needed or not and especially which version of the VMD 
model should be applied will be at least very difficult for these decay modes due to the 
lack of precise data. 

4.3 P -> /+/- 

The last of the decays via the triangle anomaly that we want to discuss is the one of a 
pseudoscalar meson into a lepton pair P — > l + l~ . As mentioned in the Chapter 3.2 a lot 
of work was done there before. 

As discussed in Chapter 3 we can give a value for the subtraction constant .4.(0), see f!3.56p . 
which contains all of the nontrivial dynamics. The results calculated with the hidden gauge 
model and the modified VMD model are compared to theoretical values and experimental 
data in Table ELT01 





hidden gauge 


modified VMD 


m 


exp. data [78] 


AM 2 


= 0) 


-21.42 ±0.04 


-21.56 ±0.06 


-21.9 ±0.3 


-18.6 ±0.9 


AM 2 


= 0) 


-5.42 ±0.02 


-5.57 ±0.06 





Table 4.10: Subtraction constant A{q 2 = 0) for decays into electrons and muons (denoted by 
the index e and /i). 

In [12] the authors presented the values of .4.(0) and the resulting branching ratios calcu- 
lated for various phenomenological models. All these results basically fall into the same 
region (the predictions from the quark models are even higher in absolute magnitude, vary- 
ing between —22 and —24.5), higher in magnitude than the experimental data. Our result 
matches very nicely the theoretical predictions, but does not overlap with the experimental 
data. 

7T° -)• e + e~ 

This situation persists for the branching ratios of the decay 7r° — > e + e~ . We present the 
unitary bound as discussed in Chapter 3 and also the branching ratios with respect to the 
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total decay rate and the one into two photons in Table 14.111 





unitary bound 


hidden gauge 


modified VMD 


T n o_ 


±e+e-l? total (lO" 8 ) 


> 4.69 


6.38 ±0.01 


6.33 ±0.02 


T n o_ 




> 4.75 


6.45 ±0.01 


6.41 ±0.02 



Table 4.11: Unitary bound and branching ratio for the decay 7r° 
different VMD models. 



e + e calculated with 



For the decay ir° — > e + e an overview of theoretical predictions calculated via different 
models and a comparison with the latest KTeV result [78] can be found in [68], see Figure 
EDI 
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Figure 4.10: Theoretical values calculated with VMD models [TjH ESI ED], with ChPT [HR E3 
[83] . modern calculations by [TJl [T2] (CLEO), and recent experimental data [75] of the KTeV 
Collaboration of the branching ratio BR(ir° — > e + e~) according to |68j . 

Most of these theoretical values can not describe the data ([78]) and this is the same 
for our calculations. The values calculated via VMD models ([73 ESI ED]) are basically 
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consistent with each other and also with the values of modern calculations [T2"l [T4"] . The 
overview contains calculations of ChPT [SU IS21 E2]- The calculated values of [SI] have 
very large uncertainties and are therefore the only values, which can describe the data. 
The values given in [82] and [83] have smaller uncertainties and an overlap with the values 
of [81]. Though these values do not overlap with each other. The authors of [82] provided 
a consistent description of the decays n° — > e + e~, r\ — > fi + fi~, and K L — > Accord- 
ing to [S3] these calculations did not take contributions of 1/N C suppressed counterterms 
into account. In contrast the calculations of [83] contain these contributions. Note that 
the values of [83] are consistent with the values of VMD models and of modern calculations. 

We will compare our results to the recent calculations given in [T2] and [2] and the 
existing experimental data explicitly, and also show briefly the values of the various ChPT 
calculations. 

The theoretical calculations of [12] and [H] and experimental data are given in Table 14.121 





m 


m 


exp. data [?S] 


Fn°-+e+e- A "total ( 10 & ) 


6.23 ±0.09 


6.26 


7.49 ± 0.38 



Table 4.12: Theoretical values and experimental data of the branching ratios of the decay 
7T —7- e + e~ . 

Our calculated values agree with these theoretical predictions. But although our value 
is a little bit higher, the theoretical predictions are still approximately three standard 
deviations lower than the experimental data. 
The branching ratio calculated by [SI| is the following: 

T^ e+e -/T total = (7± 1) x 10" 8 . (4.2) 

As mentioned this value overlaps with the experimental data, but has a very large uncer- 
tainty. The braching ratio of [S2] is 

T^ e+e -/T total = (8.3 ± 0.4) x lO" 8 . (4.3) 

This value is higher than all other calculated values and also higher than the data. The 
values of [S3] are given with respect to the decay into two photons and are consistent with 
our values: 

i\o^ e+e -/r>^ 77 = (6.2 ±0.3) x 10^ 8 . (4.4) 

t] -> i+r 

We will now discuss the situation in the //-sector. In addition to the decay r] — > e + e~ also 
the decay r\ — > fi + fi~ is possible. Our calculated values are given in Table 14.131 
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unitary bound 


hidden gauge 


modified VMD 


r r; ^ e +e-/r \otal (lO 9 ) 


> 1.78 


4.68 ±0.01 


4.65 ±0.01 




> 4.51 


11.89 ±0.02 


11.84 ±0.03 


r^ M+M -/r tottt i (io- 6 ) 


> 4.36 


4.87 ±0.02 


4.96 ±0.06 


r"jj — JU — / — >-yy (10 6 ) 


> 11.04 


12.40 ±0.02 


12.62 ±0.95 



Table 4.13: Unitary bound and branching ratios for the decays r\ — > e + e and rj — > /j> + fi 
calculated with different VMD models. 







m 


exp. data [331 EH EE] 




e~ total 


(io-*) 


4.60 ±0.09 


5.24 


< 2.7 x 10 4 




^ [i~ /r total 


(io-) 


5.12 ±0.276 


4.64 


5.8 ±0.8 



Table 4.14: Theoretical values and experimental data of the branching ratios of the decays 
rj — > e + e~ and rj — > /i+ fj,~ . 



We can compare our values to the theoretical ones and the experimental data, see Table 

EH 

Our results fall between the two different theoretical calculations given in [12] and [14"] . 
It is interesting to see that the approach given in [14], which is a little bit closer to the 
experimental data of the decay tt° — > e + e~ gives now worse predictions for the decay 77 — > 
fi + H~ if we compare it to the experimental data. As in [12], we also reach a small overlap 
with the experimental data. For the decay 77 — > e + e~ all values meet the experimental 
bound. 

The authors of [81] gave predictions for the decay 77 — > e + e~: 

^e+e-A "total = (5 ± 1) X 10~ 9 . (4.5) 

One can see again the large theoretical uncertainty. The authors of Ref. [S2] calculated 
this branching ratio as follows: 

T^ e+e -/T total = (5.8 ± 0.2) x 10~ 9 . (4.6) 

In Ref. [83] there are values given for both decays 77 — > e + e~ and 77 — > fi + fi~: 

r^ e+e -/r^ 77 = (1.15 ±0.05) x io~ 8 

rW M -AWy = (1.4 ±0.2) x 10" 5 . (4.7) 



The values of [M] are again consistent with our calculations. 
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7/ -> z+z- 

The values of the //-decays are given in Table |4"T51 Experimental data do not exist. 





un. bound 


hidden gauge 


modified VMD 


H21 


m 


r^-> e + e -/r tot (10- 10 ) 


> 0.36 


1.154 ±0.058 


1.147 ±0.063 


1.178 ±0.014 


1.86 


r»y'— s>e+e^ /r?7'— >77 (iO -9 ) 


> 1.72 


54.85 ±0.29 


54.19 ±0.42 






r^^+ M -/r to t (10- 7 ) 


> 1.35 


1.17 ±0.07 


1.14 ±0.13 


1.364 ±0.010 


1.30 




> 6.38 


5.45 ±0.15 


5.34 ±0.11 







Table 4.15: Unitary bound and branching ratios for the decays rf — > e + e and rf — > jX^ fi 
calculated with different VMD models. 



It is remarkable that the branching ratio for the decay rf — > fi + fi~ is lower than the lower 
limit given by the unitary bound as can be seen in our calculations and the ones of |14j . 
The reason is the additional imaginary part which occurs because one of the intermediate 
vector mesons can go on-shell. For the 7/ decays our calculated values are lower than both 
of the other theoretical values. 

In the r/' sector there are only ChPT calculations done by [52]: 

r v ^e+e-/T total = (1.5 ± 0.1) X 10" 10 

r^+^/Ttotai = (2.1 ±0.3) x 10" 7 . (4.8) 

The branching ratio for the decay rj — > e + e~ is very close to other theoretical values, while 
the one of 7/ — » fi + fi~ is higher. 

Summary 

In the 7r°-sector our calculated values are a little closer to the experimental data, but still 
three standard deviations lower. Our values differ slightly from the ones of [12] and [H] 
in the 7/- and 77' sector, but are still in their range. A preference of any model can not be 
given. 

4.4 P — > 7T + 7r~7 and P — > n + 7Y~e + e~ 

In this Section we will present the final results for the branching ratios of the decays 
P — > 7r + Tc^-f and P — > TT +, K~e + e~ for the hidden gauge and the modified VMD model. 
There exist several theoretical calculations ([22], [19], [21]) and experimental data ( |34] . 
[T7] . [T6]). especially in the rj sector, to compare with. We will also give predictions for the 
CP violating terms and the electric terms. Finally an upper limit for the model coefficient 
G will be given. 
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77 — > 7T + 7r 7 and 77 — > 7T + 7r e + e 

We start with the dependence of the rate of the decay rj — > 7r + 7r~e + e~ on the invariant 
mass of the decaying particle. We plot the leading (magnetic) term calculated, with the 
hidden gauge model and the modified VMD model, and the electric term which is scaled 
up to make it comparable (see Figures 14.111 and I4.12p . 




Figure 4.11: Dependency of the differential decay rate of the decay rj — > 7r + 7r e + e on the 
change of the invariant mass of the pions. 




The behavior is quite similar to the graphs we showed before. With respect to the depen- 
dence on the invariant mass of the electrons, there is a large peak at low energies. The 
behavior of the invariant mass dependence of the pions is much broader with a smaller 
peak. 

For the decay 77 — > 7r + 7r~/i + /j _ , the graphs are shown in Figures [4.131 and 14. 14[ 

It is interesting that in these decays the electric terms seem to be much bigger than in the 

previous case. But this is actually not the case if one compares the calculated values of 
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the decay rates and branching ratios. 

These values for the r/ sector are listed in Table 14.161 Here the branching ratios relative 
to the total r] decay width and to the width of the r\ — > 7r + 7r~7-decay are listed. The 
asymmetry term is defined according to Eq. (2) of [23] with the mixed term normalized 
to the total rj — > TT + ir~fi + ^~ width, see Eq. (I3.90P of Chapter 3. The electric terms are 
independent of the various VMD models and normalized to the total width. 
One can see that the branching ratio of rj — > 7T + ir~ fi + fi~ is much smaller than the other 
branching ratios. This is due to the smaller phase space. The electric terms are about 
two orders of magnitude smaller for the r] —> n + 7T~e + e~ decay and one order of magnitude 
smaller for the rj — > 7i + tt~ fi + fi~ decay, even if we set G to one. So it is obvious that the 
magnetic terms still determines the leading contribution. The difference in predictions 
of the various VMD models are highly visible for the case of rj — > 7r + 7r~7 and the rj — > 
7T + 7r _ e + e~-decay. It is interesting that the CP- violating asymmetry term has a negative 
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hidden gauge 


modified VMD 


7] — > 7T + 7T 7 


T' 7T + 7T j/^total 


(10- 


- 2 ) 


4.97 ±0.27 


4.79 ±0.19 


7] — > 7r + 7r e + e 


~^Tr+ir~ e+e - /h total 


(10- 


- 4 ) 


3.14 ±0.17 


3.02 ±0.12 


~^Tr+Tr~ e+e - /T tt+tt— 7 


(10- 


- 3 ) 


6.32 ±0.01 


6.33 ±0.01 




(10- 


- 2 ) 


-3.88 ±0.11 


-3.95 ±0.08 


E — terms 


(10- 


- 6 ) 


1.6 G 2 independent on VMD 


7] —7- 7T + 7T fl + fl 


h tt+tt ^t+/i /^total 


(10- 


- 9 ) 


8.65 ±0.39 


8.64 ±0.25 




(10- 


- 7 ) 


1.74 ±0.02 


1.81 ±0.02 




(10- 


- 2 ) 


-1.28 ±0.03 


-1.28 ±0.02 


E — terms 


(10- 


9 ) 


0.889 G 2 independent on VMD 



Table 4.16: Branching ratios for the decay r] — > tt + tt 7, r] — > tt + tt e + e and r] — > tt + tt /i + ji 
calculated with different VMD models. 



sign. This is due to the negative sign of the VMD term relative to the contact term. In the 
case of the magnetic term this term appears squared, so it does not affect the sign there. 
We can now compare our values with other data, see Table 14.161 and 14.171 Note that the 
authors of [19] and [21] used the hidden gauge model as well, while [22] calculated in the 
framework of unitary chiral perturbation theory. So it will be more interesting to compare 
our values with [22], because an agreement would be a verification of both models. The 
discrepancies between our values calculated via the hidden gauge model, and the ones given 
in [IH] for the r] — > 7i + n~'~f and in [21] for the rj — > ir + 7r~e + e~ decay can be explained by 
the modern values of the constants fo and as well as the 77/7/ mixing angle. The 
data of [17] are the most recent published data and also the most interesting. 

For the rj — > 7r + 7r~7-decay all theoretical values represent the data very well. The modified 
VMD result is closer to [22] than the hidden gauge result. 

The decay r/ — > n + 'K~e + e~ is more interesting. One can see by comparing Table 147161 and 
14.171 that the hidden gauge model does not agree with [22], but that the modified model 
gives nearly the same value. The experimental situation is even more interesting. The 
values of the hidden gauge model actually represent the data of [31] better than the modi- 
fied model, although both are consistent. The recent data of the KLOE measurement [17] 
are much smaller than the other experimental data and also have extremely small errors. 
Neither other theoretical calculations [22] nor our values have an overlap. It would be very 
interesting to see whether the measurement of the WASA@COSY experiment is compatible 
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ESI, EI] 


PDG 01] 


KLOE P2] 


?7 —7- 7r + 7r 7 


r/r tota; (io- 2 ) 


4 - b8 -0.09 


5.22 


4.6 ±0.16 




7] —7- 7r + 7r e + e 


r/r tota , (lo- 4 ) 


9 qq+o.08 

Z - Jy -0.11 


3.2 ±0.3 


4.2 ± 1.2 


2.68 ±0.09 ±0.07 


r/r 7r+7r - 7 (io- 3 ) 


6-39+^1 




9.2 ±2.5 




^ (io- 2 ) 






-0.6 ±2.5 ± 1.8 


7] —7- 7T + 7T 


r/r iotai (10-9) 


7 c;+ 4 - 5 

' ' d -2.7 




< 3.6- 10 5 




r/r 7r+7r - 7 (io- 7 ) 


i «i +0.95 
i - ui -0.55 







Table 4.17: Theoretical values ( [22], [E] and [21]) and experimental data ([M] and [T7J) of the 

branching ratios of the decay r\ — > -K + iT~~f, r\ — > ir + -K~e + e~ and r] — > ir + ir~ fj, + fi~ . 



with the one of KLOE, when their new data appear. (The most recent published data of 
WASA@CELSIUS are already included in the PDG value and can be found in Ref. |16j.) 
If the data were confirmed a new theoretical view on this decay channel would be needed. 
The ratio r 7r + 7r - e + e -/r 7r + 7r - 7 is relatively constant within the various VMD models, which 
makes sense, since the VMD factors should approximately cancel. The errors are smaller 
because some prefactors cancel also. Note that our values are in agreement with the ones 
of [22]. 

The comparison of this ratio to the recent KLOE data is again of special interest (see Table 

an. 



(10-3) 



KLOE 
CLEO 



6.77 ±0.44 



KLOE 
PDG 



5.83 ±0.31 



PDG 



9.2 ±2.5 



Table 4.18: The branching ratio T w + w - e + e -_/T w + n -^ calculated with different data of the KLOE 
and CLEO experiment ([IT] and (SO) an d the PDG value ([53]). 



Normalized to the PDG value the fraction is smaller than the theoretical predictions and 
has no overlap with these values. The normalization to the recent CLEO data ([85]) gives 
a larger value, because the data of the CLEO experiment is smaller than the PDG value 
(BR(r] -> 7T+7T-7) = (3.96 ± 0.14 ± 0.14) x 10~ 2 ). This ratio has a small overlap with 
the UChPT value ([22]), as already pointed out in Ref. [17], but barely overlaps with our 
results. Note that the PDG value does not include the CLEO data. The PDG value is larger 
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than all the other values. It has no overlap with any theoretical value. A measurement of 
this branching ratio in one experiment would be usefull. 

We can now take a look at the CP violating asymmetry term A^, see Eq. f)3.90p of Chapter 
3. Theoretically it was estimated by Ref. (23] as A^ = 2.0- 10~ 2 G which is about half of our 
calculated value. A^ was measured for the first time via the forward-backward asymmetry 
by Ref. [17]. According to these data we can give constraints for the model coefficient G. 
To make sure that our calculated values are consistent with the KLOE data is has to be 



-0.9<G<1.2. (4.9) 
This agrees with the assumption that G < 0(1) is natural. 

For the T) — > 7r + 7r~/! + /i~-decay, the choice of the VMD model does not have a serious effect 
on the branching ratio, which is remarkable. Our values are consistent with the ones of 
[22] and meet the upper experimental bound easily. 



rf — > 7T + 7r 7 and rf — > 7r + 7r e + e 

In the curves of the differential decay rate of the decay 7/ — > 7T + 7T~e + e~ there is a very 
strong contribution of the width as one can see in Figure 14.151 For the leading mag- 
netic contribution the curve is shifted to the right for larger vector meson masses. The 
dependence of the asymmetry term is very interesting. 




For the branching ratios the situation is somewhat different as we will discuss using the 
calculated values of the branching rations, see Table 14.191 

The calculation of the modified model now gives larger values than the calculations of the 
hidden gauge model. We can see that the CP violating terms are one order of magnitude 
lower than we the ones in the 77-decays and are of positive sign. The electric terms are 
visibly smaller than in the 77-decays. As we pointed out before we can assume G < 0(1) 
and so the electric terms do not have any effect. 
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hidden gauge 


modified VMD 


th. values [22] 


exp. data [M] 


77' —7- 7T + 7T 7 




0.294 ± 0.027 


0.308 ±0.016 


n 294+ - 027 




rf —> 7T + 7r e + e 


r7r+^-e+e-/rtota« (lO 3 ) 


2.17 ±0.21 


2.27 ±0.13 


9 1 O+0.19 
z - io -0.32 


9 r+1.3 

z - d -1.0 


r7r+7r _ e+e _ /r-Tr+Ti-— (lO ) 


7.368 ± 0.043 


7.370 ± 0.059 


7 9 4 +o.09 

' - z ^-0.15 




^/G (10- 3 ) 


2.83 ±0.29 


2.02 ±0.18 




E - terms ( 10_7 ) 


1.67G 2 


0.85G 2 




77' —7- 7T + 7T fl + fi 


^ -K+-K- ii+ n~ total (10 5 ) 


2.20 ±0.30 


2.41 ±0.25 


L -° ' -0.75 


< 23 


r7r + 7r _ /i+/l _ /Ttt+tj- — -y (lO ) 


7.46 ± 0.32 


7.82 ± 0.42 


5 4+3.6 




(10- 3 ) 


4.41 ±0.33 


3.03 ±0.19 




E - terms (l°~ 8 ) 


1.28G 2 


0.65G 2 





Table 4.19: Branching ratios for the decay rj' — > tt + tt 7, 7/ — >• 7r + 7r e + e and rj' —> tt + tt /i + ji 
calculated with different VMD models, other theoretical values and experimental data. 



In contrast to the ?7-decay, the agreement of our values and the ones of [22] does not hold 
any longer. For the 77' — > 7r + 7r~7 and the 77' — > 7r + 7r~e + e~ channel the discrepancies are 
very small. However, for the decay 7/ — > n + 7i^ fi + fi~ there is only a small overlap, although 
the errors of [22] are very large. 

Both of our VMD models as well as the values of [22] are consistent with the data [34], 
though our values are less constraint by the data. More precise measurements, especially 
for the 77' — > 7r + 7r~ fi + fi~ decay, are needed to see which model represents the data the best. 



Summary 

For the 77 decays the values calculated with the modified VMD model are closer to the 
experimental data. One can see that the modified model is indeed useful in the box anomaly 
sector. The lesser agreement for the 77'-decays is due to the scarce data. This could also 
be the reason why the values of [22] are different especially for the decay 77' — > 7r + 7r~ fi + fi~ . 
The asymmetry term was only measured for the decay 77 — > 7r + 7r~e + e~ . Although the 
errors given for the total branching ratios are extremely small, the ones for the asymmetry 
term are still much higher than the value itself. For the other decays we made predictions 
for the asymmetry term, which will be more accurate if the parameter G can be given more 
precisely. 
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Chapter 5 

Summary and outlook 



In this work we studied anomalous decays of pseudoscalar mesons. We calculated explicit 
expressions for the decay rates and branching ratios and discussed the relevance of the 
form factors of the various vector meson dominance models. Thereafter, we presented 
our results of the branching ratios for the different vector meson dominance models and 
compared them to other theoretical values and experimental data. 

For the decay P — > Z + Z~7 our results represented the experimental data very well. Also 
the theoretical values calculated by other groups totally agreed with ours. The difference 
between the values calculated with and without VMD models was very pronounced for the 
decays 77/77' — > fi + fi~'~f. We were able to conclude that the VMD factor is needed here to 
represent the data, but we had no preference on one of the models. For the rj' — > Z + Z~7 we 
also showed the contribution of the width. 

The situation in the P — > / + /~/ + /~-channel was slightly different. Our values and the other 
presented theoretical values of the decay tt° — > e + e~e + e~ could represent the experimental 
data equally well. For the interference term we found differences to the results of the older 
calculations of [7] which had also been found by [8] . We calculated very different values for 
the decays 77 — > l + l~l + l~ in comparison to the other theoretical groups. We can refer again 
to [S] where an error in the calculations of [7J was found. Unfortunately experimental data 
are very scarce and only give upper bounds that all theoretical predictions could meet, 
so we could not constrain our results there. Theoretical calculations and experimental 
data for the respective 7/-decays do not exist so that we could only make predictions. In 
the 77' sector we also found contributions of the width in the decays 77' —> e + e~e + e~ and 
77' —> /i + /i~e + e~. 

Our improvements in the case of P — > were very limited. We basically followed the 
calculations of [12] and [H]. Our calculated values for the decay n° — > e + e~ are very 
close to the ones they calculated and still far away from the experimental data. For the 77 
decays our results fall between the values of [12] and [2] and even have an overlap with 
the experimental data for the decay channel. The values for the 77' decays are again 

very close to the ones of [12] and [H], but there exist no data to compare with. 

In general, we found that the differences between the VMD models in the triangle anomaly 
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sector are insignificant. This result is welcome, because the hidden gauge model could 
already describe the data very well. 

In the box anomaly sector, we found that the modified VMD model led to improvements. 
The calculated values for the branching ratios of the decay 77 — > ir + iT~^f represented the data 
very well for both VMD models, but the improved VMD model achieved a better agreement 
with the theoretical values calculated via unitary chiral perturbation theory ([22])- This 
was also valid for the decay 77 — > n + 7r~e + e~ . On the other hand the experimental situation 
is different, because the recent measurements done by KLOE supplied distinctly smaller 
values. Here the modified VMD model was closer to the experimental data and could 
almost reach an overlap. The calculated values for the decay 77 — > tt + tt~ fi + fi~ matched 
with the ones of other theoretical calculations and met the upper experimental bound. 

In the decay 77 — > 7r + 7r _ e + e~ we also analyzed the CP-violating asymmetry term and 
compared our values to recent experimental data. We calculated an upper bound for the 
model specific factor G. This factor scales a CP-violating, flavor-conserving local four- 
quark operator which is sensitive to the ss content of rj and 7/ and determines the strength 
of the additional electric form factor E. We verified the claim that it could be of natural 
size. We also gave a prediction for the asymmetry term of the decay 77 — > tt + tt^ fi + fi~ . 

For the 77' decays the situation was again different. The only existing data were the ones 
of the decay 77' — > ir + iT~e + e~ . Here all models represented the data very well. For the 
decay 77' — > 7r + 7r~/i + /i~ our values were very different from the one calculated by [22], but 
the experimental situation is too scarce that a preference of any work can be justified. We 
found again very interesting contributions of the width in the 77' sector. Especially for the 
asymmetry term we had a change in the algebraic sign in the region of the vector meson 
mass. 

We completed the calculations concerning the anomalous decays of the 77 and the 7r° in 
the framework of vector meson dominance. The only decays where the experimental data 
could not be described by theoretical calculations were the ones into two leptons P — > l + l~ 
and the decay rj — > n + 7T~e + e~. If the experimental data are confirmed, the probability 
that new physics is needed to describe the data is high indeed. The same holds for the 
decay rj — > 7T + n~e + e~, although the theoretical values and data are much closer than in 
the decay into two leptons. 

The anomalous 7/-decay 77' — > 7r + 7r~7r + 7r~ should be investigated. The kinematics should 
be very similar to the ones of the decays P — > especially an interference term 

will contribute, but the form factor would be of special interest. 



Appendix A 
Kinematics 



A.l The parallel boosts 

Throughout we assume that all boosts are performed parallel (or anti-parallel) to the 
virtual- or real-photon axis which should point parallel to the z-axis. Therefore only the 
th and 3 rd components of the 4- vectors will be affected by the boosts (e.g. q° and q z ), 
whereas the 1 st and 2 nd components q x and q v will remain untouched. This is the reason 
for the notation qn and q^ introduced in (I2.49P and (I2.50p . Without loss of generality, we 
can therefore always mangage to rewrite (13.611) as 








( P°- \ 








H 




V p-h J 






( p°~ \ 




H 


-p+± 






V p-„ J 






(a.i; 



(A.2) 

Note that the frames and boosts have been chosen in such a way that Pj_ = p_i_ = = 0j_. 
Therefore, we always have 

P-i = -P+± , ( A -3) 
k +± = -k_ ± (A.4) 

as it was applied in (IA.2j) . 



A.2 The relevant frames 

The coordinates and four- momenta in the P rest frame are denoted here by a tilde (~), the 
ones in the p+P- rest frame by an asterix(*), and the ones in the k + k_ rest frame by a a 
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diamond (°). 

Thus, the relation flA.ip is given in the P rest frame by 





ft \ 




( K > 






p+± 


H 


-p+± 


H 


Ox 


?n ) 




v p--, ) 







{ ~ k ll J 



whereas in the rest frame it reads as 






j pp 



0± 

0|| 





f P%° \ 




< * \ 


H 




l + l 






V -p+*„ J 








k* 





(A.5) 



(A.6) 



Switching the asterix(*) with a diamond (°) and p with fc this relation holds in the k + k^ 
rest frame. 



A. 3 Comparison with other kinematics 

Note that 9* + is exactly the angle 9 defined below Eq. (6) of Ref. [72]. Moreover, since 
the three-momentum P* of the P meson in the p+p- rest frame is identical to the three- 
momentum k* in that frame (see (I A - 6 [) ) . the angle 9* + is identical to the angle 9 V defined 
below Eq. (1) of Ref. [22], i-e. identical to the angle between the p + three-momentum and 
the P three-momentum in the p+P- rest frame: 



9* p+ = 9 of Ref. 



9 n of Ref. [23]. 
Finally note 

(p +± ) 2 = (p + ) 2 sin 2 ^ + = ( P ;) 2 sin 2 ^ + = ( P ;) 2 sin 2 ^ 
and in total analogy 

o 

i 2 „;„2 no — /i,o \ 2 • 2 



(k- 



k_ sin 2 ^- = (k!fBin a 0£_ = (1st) sm 2 9 k . 



(A.7) 



(A.8) 



(A.9) 
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See Ref. [23] for the definition of 9^: i.e. df. is the angle between the fc_ three-momentum 
and the P three-momentum in the k_k + rest frame. Note that the P three-momentum in 
the k-k + rest frame points into the negative ^-direction, therefore 9k should be replaced 
by 7r — 9k if a common coordinate system (which involves 9 P ) is used. Equations (1A.7j) . 
( ]A.8j) and flA.9j) are the essential formulae. 



The kinematic of the decay P(P) — > pf(p + ) P2 (p~) p% (k + ) p~l (&_.) is very close to the kine- 
matics of the K14 decay, the decay process of one kaon into two pions, one anti-lepton and 
one (to the lepton corresponding) neutrino, e.g.:K + {p K ) — > ix + (p + ) n~ l + (ki) i>i(k v ) . 
The corresponding kinematics were introduced and described in Cabibbo and Maksy- 
movicz [51], Pais and Treiman [S7J and Bijnens et al. [7T]. There is, however, one subtlety: 
the mass of the neutrino v% is of course assumed to be zero, whereas here both leptons have 
the non- vanishing mass m e . This induces changes in e.g. (lA.20j) . (1A.21j) (IA.25j) . (IA.26j) . 
( IA.27P see below. 

Refs. [SS], [HZ] and [7T] utilize the follwing five variables (see Fig. 1 of Ref. (HS| or Fig. 5.1 
of Ref. [7Tj ) transcribed to the decay P — > p\P2PzPi- 




Figure A.l: Kinematics of the Kn decay [7T] 

1. s pp = (p + + p-) 2 (corresponding to R 2 in [86] and s n in [871 E]), the effective mass 
squared of the P1P2 system, 

2. Skk = (k + + k-) 2 (corresponding to K 2 in [86] and si in [HJJ EI]), the effective mass 
squared of the p 3 p 4 system, 

3. 9 P (corresponding to 9 in [SB] and 9 n in [H7J EE]), the polar angle (O^^p^vr) of the 
Pi{p + ) i n the P1P2 rest frame with respect to the direction of flight of the P1P2 in the 
P rest frame, 

4. (corresponding to ( in [86] and 61 in Refs. [UJI7T]), the polar angle (0<9k<ii) of 
the Pi(k~) in the P3P4 rest frame with respect to the direction of flight of the p^Pi in 
the P rest frame, and 



s2 



APPENDIX A. KINEMATICS 



5. tp (<fr in [SU [7T] and tp in [ST]), the azimuthal angle (0 < <p < 27r) between the plane 
formed by in the P rest frame and the corresponding plane formed by the P3P4. 

Ref. [87] uses a convention for the metric opposite to the standard Bjorken-Drell metric. 
The latter is applied here and in Refs. [86j [71] . Thus 



9 P = 9oi Ref. [7J = 9 n of Ref . [23] = vr - of Ref . jHSl 
= 7T - ff of Ref. [STJ d] = % - 12 of Ref. [S| 



and 



^ fc = ^ e of Ref. [23] = C of Ref. [86] = ^ of Ref. [87J [JT] . = fl 34 of Ref. [87J [JT]. 



In Ref. [861I7T] the following explicit construction can be found: Let v be a unit vector along 
the direction of flight of the 7r + 7r~ in the P rest frame (Sp), i.e. v = — k. Furthermore, 
let c (d) a unit vector along the projection of p + (k + ) perpendicular to v (— v), 



c 
d 



Then one has 

cos 9 p 
cos 9k 
cos<p 
sin ip 



P+ ~ v(v P 
k+ - 



f (v • p 4 
vfv • kj 



VK + ( v ■ ^ 



= k+, 



-k_ 



(A.10) 
(All) 



—v • p + = k ■ p + = cos 9* , 
—v ■ k + = k ■ k + = — k ■ k = 
c • d = p +± • k +± = -p +± • k_ 
(c A v) ■ d = -v ■ (c A d) 

A f— k_, ) I = -sin<^ 



k- p 



'+± 



cos 91 , 



cos< 



(A.12) 

(A.13) 
(A.14) 

(A.15) 



Thus ip = 7T + (f). 
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A. 4 Invariant expressions of the decay momenta 

The interesting relations read: 





p 2 = 


ip++P-] 


f 


= 2m 2 + 2p + ■ p^ 


- S pp j 


(A.16) 


(p+ 


+P-) 


■ (P+ - P- 


) 


= m 2 — nip = 




(A.17) 






(P+-P-. 


) 2 


= 2m 2 , - 2p + ■ p_ 


- Am p Spp Spp(3p , 


(A.18) 




k 2 = 


(k+ + k_] 


) 2 


= 2m 2 k + 2k + -k. 


- s kk ) 


(A.19) 


(*+ 


+ k_) 


■ (k+ - k_ 


) 


= m 2 k -m 2 k = 0, 




(A.20) 






(k+-k-] 


) 2 


= 2m\ - 2k + ■ k. 


_ = Am\ ~ s kk = -Skkfit , 


(A.21) 


0+ 


+ p- H 


- k + + kJ 


) 2 


= {p + k) 2 = P* 


= m 2 P 


(A.22) 


(p+ 


+ P-) 


■ (k + + k. 


-) 


= P-(k + + k_)- 


- s kk = (p + +p~) ■ P- s pp 












= 2 ( m P Spp ~ 


Skk) 


(A.23) 


(p+ 


-P-) 


■ (k + + k. 


-) 


= (P+~P~)-P 














Pp 2^ ^ ^pp, 


m 2 P ,s kk ) cos 9 P , 


(A.24) 


(p+ 


+ P-) 


■ (k+ - k- 


-) 


= P-(k+-k-) 














= +/3 k 1 -X 1 / 2 (s kk , 


m 2 p,s P p) cos9 k , 


(A.25) 



(p + -p_) • (k + -k-) = -±(3p(3 k (m 2 P - Spp- s kk ) cos 9 P cos 9 k 

+ y/s p ~p'f3p^/s k ~ k 'f3 k sin 9 P sin 9 k cos (f> 

(A.26) 



e^{p + -P-Yp u {k + -k.) a k^ = -mpe^ k k\p + -p^{k + -k^) k 



g-^ ^ (jripi Spp, s kk )-^y Spp(5p\J s kk (5 k 
x sin^p sm9 k sin0 

[\(m 2 p,s P p,s kk )\(spp,m 2 p ,m 2 p )\(s kk ,m 2 ,m 2 )] 1/2 . Q . Q . , Q „s 
— — — — sin Up smu k sin (p. yP^..Z(j 

V ' SppSkk 

Some of the relations given above deserve some additional remarks: Equation (1A.23|) can 
be rewritten as 

k ■ (p+ +P-\ = |( m p - Spp - s kk ) . (A. 28) 

It can be combined with correspondingly rewritten equation (IA.24|) 

fc-(p+-J>-) = -k* (p* -pi) =-2|k*||p;| cosf), 



V Spp 

Spp, m P , s kk ) cos9 p , (A. 29) 
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such that the following relations hold: 



k-p+ 
k ■ p_ 
p-k + 
p ■ fc_ 



4 ( m p s pp s kk) 4/3pA ^ {spp,Trip, Skk) COs9 p , 
\ (mp - Spp - s kk ) + \P P \ 1/2 (s P p, rrip, s kk ) cos9 p , 

°pp 

Spp, m%, s kk ) cos6 k , 
Spp - s kk ) - \p k \ l/2 (sp P , nip, s kk ) cos9 k . 



IK 



(A.30) 
(A.31) 
(A.32) 
(A.33) 



Note that these expressions deviate from Eqs.(lO) and (11) of Ref. [23], where there is a 
factor =Fl/2 instead of Tl/4 in front of (3 P . In fact, the normalization of Eq. flA.29j) exactly 
agrees with the normalization of the first relation of Eq. (3) from Ref. [57], when Eqs. (2') 
and (3') of that reference are inserted into this relation. 
Using (1A.25|) and (1A.26|) we will get the following expressions: 



(k+ - k_) -p + 



(k + - - p. 



(P+-P-) ■ k+ 



(P+-P-) ■ k- 



\ 1/2 (s pp , rripp 2 , s kk ) - (m 2 P - s pp - s kk ) f3 p cos0 p 
+ \ y/Sp~p~f3 P y/skk(3k sin 9 P sin 9 k cos , 



X 1/2 {s pp , mp, s kk ) + (m 2 P - Spp - s kk ) (5 P cosft 

- ^y/^PpVskkPk sin 9 p sin 9 k cos </> , 

~\ ^ 1/2 ( S PP 5 m Pi s kk) + K ~~ S PP ~ S kk) Pk COs9 k 

+ ly/^Ppy/skkPk sin 9 P sin 9 k cos <p , 
~\ \ l/2 (s P p,m 2 p,s kk ) - (m 

- \^/~s^l3 P ^/~Skkfik sin 9 p sin 9 k cos 



Spp - s kk ) (5 k cos9 k 



■ (3 k cos 9 k 
k cos 9 k 

■ (3 P cos 9 p 

■ (3 P cos 9 P 

(A.34) 



In order to derive the relation (I A. 26[) we used the relation 

|k*||| = — (x 1/2 {s kk , m 2 P , Spp) - (m 2 P - s pp - s kk )(3 k cos9 k j . (A.35) 

4 v s pp V / 



This expression can be calculated via 

|k*|| = |kl| 2 - \k_\ 2 ± with |kl| = y / (^ ) 2 

The zero component /c*° can be expressed via (IA.23|) and (1A.25|) to 

1 



nr. 



kf 



4,/F 



m. 



vv 



s kk ) ± f3 k \ 1/2 {s kk ,m 



pi s pp) 



cos 9i 



(A.36) 
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With this relations we can calculate the expression 



(p* + -p*_ )„ ■ (n-ki)n -( P+ -p_) ± ■ (k + -k_ 



2p 



+11 



k^-2k* 



-2 ( P ; • k) (k* • k) + 4(p* • k)(ki ■ k) + 4p^k_^ 

Eg. JA~T2l Eg. fA~l41 

— 2|p* | cos# p |k*| +4|p* |cos0 p |k* ||| 
+ 4|p+| |k* | sin 6 P sin 9 k cos 



(A.37) 



which is essential for the structure of ()A.26j) . 
Finally we list the scalar products: 



k+-p- = v^ 1/2 (s pp , m 2 P , s kk ) [P k cos6 k + /3 P cos6 p ] 

o 

+ o { m P ~ S PP ~ Skk ) t 1 + fikPp COS ® k COs6 W 
o 

- \y/s^P P \fskkPk sin 9 P sin k cos , 
k--p+ = -]-\ 1/2 (s P p, m P , s kk ) [f3 k cos 9 k + (3 P cos9 p ] 

o 

+ I i m P ~ s pp ~ Skk ) I 1 + PhPp cos ° k cos ®p\ 
o 

- \y/Sp~p\P P \fskkPk sin 9 P sin k cos , 

k+ ■ p+ = ^A 1/2 (s pp , m 2 P , s kk ) [(3 k cos9 k - (3 P cos9 p ] 

o 

+ rj { m P ~ s pp ~ s kk) [1 - PkPp cos 9 k cos9p\ 

o 

+ \ y/Sppfip y/skkPk sin 9 P sin 9 k cos , 
k- ■ p- = -^A 1/2 (spp, mp, s kk ) [f3 k cos 9 k - (3 P cos9 p ] 

o 

+ rj ( m P _ S PP ~ S kk) [1 - PkPp COS fc COS^p] 
o 

+ \y/Sp~p^Ppy/sk~kfik sin P sin 9 k cos . 

(A.38) 
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A. 5 Projection tensor 



The projection tensor, which we needed in Chapter 3, can be constructed via the summation 
over the final spins s_ and s + of the current s_; k + , s + ) 

1/2 1/2 

O j _ l f 1 '{k_,k + ) = jti(k_, s_; k + , s + ) jj l ,(k_, s_; k + , s + ) 

s_=-l/2 s+=-l/2 
1/2 1/2 

= e2 Tr i u ( k -i s -)' a (k-,S-.)"f IM v(k + ,s + )v(k + ,s + )"f IM/ } 

S-=-l/2 s + =-l/2 

= e 2 Tr[(#L + m)7 M (^ + -m)7 M /] 

= 4e 2 [(A;_) AI (A; + ), i / + (k+)^)^ - g^{k- ■ k+ + m 2 )] 



-2e 2 



g^'(k + + k_) 2 - (fc + +fc_) jU (fc + +fc_) /Lt / + (k + -k_)^(k + -k_) t 



e 2 k 2 x 2 



k^k^/^ {k + k )^{k + k 

^2~ 



A- 



(A.39) 



A. 6 Decay rate 

We will now discuss the decay rate for a decay of a particle P with momentum pp into 
four particles with momenta fc + and k_. The notation of the particles follows the 

ones of the momenta. The decay rates into two and three particles are also given, (1A.43|) 
and CQ4j) . 

The decay rate is defined as follows: 

dr = | A| 2 d$ 4 (pp;p+,P-, (A.40) 

with the four body phase space: 
d$ 4 (j0p;p + ,p_,/c + ,/c_) 



4 x , . d 3 p+ dV d 3 A; + d 3 k. 

o {pp-{p++p- + k + + k^)) 



(2n)^2E p+ (2Ti)nE p _ (2vr)32£ fc+ (2n)nE k _ 



(A.4i; 



\A\ 2 is the matrix element squared. 

We investigate the case where P decays into three particles with momenta p+,p~ and k, 
where the third particle k, in our decay it is always a (off-shell) photon, decays into the 
remaining two particles k + , k_ via a two-body decay: 



d<$>±(p P ;p + ,p_,k + ,k_) = (27r) 3 d$ 2 (A;;p + ,p_) ■ d^ 3 (p P ;p + ,p_,k) dk 2 (A.42) 
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with the following expressions for the two body and three body decays 

1 1 |k*| 
16tt 2 (2tt) 4 v /s^" 

(27ry8^ |p + l|k|dmp+p -^^ 



d<S> 2 (k;p + ,p_) = ^_i__KLdnL, (A.43) 



d$ 3 (pp;p + ,p-,fc) = ,„ . qo |p;i|k|dm p+p _dQ* dQfc. (A.44) 



The decay rate becomes 

dr = |A| 2 * ^ ^ |k*| |p* + | |*| dm p+p _dQl + dQ; + dQ k dk 2 . (A.45) 

The angles are defined as df2 = d0dcos#; df2^ + is the solid angle of the k + in the k + k~ 

rest frame, df2* + is the angle of p + in the p + p~ rest frame; and dQ^ is the angle of k in the 
P rest frame. 

We use the relation d^/s p + p - = 2 , 1 dg p + p - in the integration. The momenta are the 
same as defined in the previous Section. 

Because the squared matrix element only contains the angles 6^, 9* and <fi we can integrate 
over the remaining angles. Therefore, we combine the azimuthal angles of VL <> k+ and Q p+ 
to the corresponding angle tp as the difference between the planes and integrate over the 
remaining independent angle in Fig. IA.1I The angles Q k stay untouched and can be inte- 
grated out. The decay rate becomes 



dr = \A\ 2 = _ -L^rlffl I pi I Ikl — =ds„ +p -dcos^dcos#:d</?dfc 2 . (A.46) 
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Appendix B 

Translation formulae for the mixed 
and interference term of the decay 

p -> i+ri+r 



Most of the expressions are invariant and independent of the different rest frames and 
therefore easy to be calculated in terms of the equations given in Appendix A: 



S14 



S23 



m, 



Pl2 s VZ ~ /^34 s 34 + A 1/2 (?Tip, S U , S 23 ) [Pl2 COS 9 12 + /3 3 4 COS 34 ] 



+ (m^ - S12 - S 34 )/3l 2 /334 COS 6>i2 COS 6*34 



2>/si2\/s34/5i2^34sin0i2 SU1 34 COS0 , (B.l) 
^12^12 - /3 34 S 34 - A 1/2 (m|, S14, s 23 ) [P12 COS 012 + /3 34 cos 34 ] 

+ (m 2 - sia - s u )(3 12 (3 34 , cos i2 cos 34 

-2 v / si 2 A/s 3 4/3i2/3 3 4sin0i 2 sin03 4 cos0 . (B.2) 



Using the above expressions, one can easily calculate X(m%, s u , s 23 ) as well as (3 U and /3 23 . 
The expressions for the angles 0i 4 , 23 and <p are more complicated, but can also be given 
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APPENDIX B. TRANSLATION FORMULAE FOR THE MIXED AND 
INTERFERENCE TERM OF THE DECAY P L + L~L + L~ 



in terms of the invariants: 



'ii 



cos 6 1 



Sl2 - S34 - \\ 1/2 {m 2 pi S12, S 34 ) [/3i2 COS 12 - /3 3 4 COS 0; 



31 J 



cos 2 



23 



/3i 4 A 1 / 2 (m 2 3 ,si2,s 3 4) 
S12 - s 3 4 + |A 1/2 (m 2 3 ,si2,s 34 )[/3i2 cos ^12 - /3 3 4 cos0 34 ] 



/3 2 3X 1/2 {mp, S12, s 34 ) 



(B.3) 
(B.4) 



sin #i4 sin #23 cos 1 



1 



2VsU-v/S23/?14/323 



S12 + s 34 - -(m| - S12 - s 34 ) [1 - /3i 2 /3 3 4 cos0i 2 cos0 23 ] 
si2/9i2\/s34/334 sin 9 12 sin 34 cos 



m p ~ 2^ m2p ~ 312 ~ Su ) + ^ 12 ^ 34 cos ° 12 cos ^ 34 ] 

Sl2/3l2 V^34/^34 Sin 012 sin 34 COS ) /9i2/934 COS 012 cos 



734 



(B.5) 



Appendix C 



Further calculations for the decay 

P -> /+/- 



C.l Numerator of the P — )> l + l amplitude Eq. (13.481) 



To calculate the numerator of the Integral in A we first deal with Tr [v (p', s + )u(p, s_)] 
in (13.481) . This form is rather unusual, because normally the completeness relations have 
a form with the same particle operator and the same momenta, e.g. u(p, s)u(p, s). The 
derivation of such terms can be found in Appendix A of [61]. Therefore we calculate 
structures with different momenta (u(k, s)u(p, s')) first and insert the different operators 
via the expressions v(k,s) = j5u(k, — s), respectively, v(k, — s) = —u(k,s)j5. A general 
form for all structures is given in equation (A6) of that paper. Using this expressions 
general forms of projection operators on the singlet and triplet state of an outgoing lepton 
pair can be calculated. Because the total angular momentum of the lepton pair is J = 
it can only be either in a singlet 1 5'o or in a triplet 3 Po state. The decaying particle is a 
pseudoscalar meson with C — +1 and P = —1, such that CP = —1. Since CP of the 
lepton pair is given by CP = (— 1) L+S (— l) s+1 = (— l) s+l the outgoing lepton pair has to 
be in a singlet state to keep CP- invariance. We can now replace the term in (I3.48P by the 
projector on the singlet state given in equation (A16) of [M]. This reads, 



P {0) (p,k) = -^=[v(p,+)u(k,-) + v(p,-)u(p,+)} 

= r (-^m7 m 75 + + P 7 s) 

2 (m p m k +p - fc) 2 

= -1= [-2m(p + fcVrV + - P p k a ) + frys)] , (C.l) 
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where 



An = mk^ + mp^, 
1 

t 



T^ u — —6 nv P<J (k p p p p k J 



andt = 2P = {p + kf = \{m 2 +p-k) (C.2) 
were inserted. With this the amplitude changes to 

A =HJ£r w= w^m, - w x VMD{k2 ' iq ~ kn (a3) 

L^y is the singlet projection of the final lepton pair, which can be given by calculating the 
trace 

V = Tr [P(°\p,p'h, + mi ) 7,] = -2iV2^e^ aT p'°q T . (C.4) 

Combining the remaining momenta and dealing with the two total antisymmetric tensors, 
we can derive 

e^VflT = = {m 2 P k 2 - (k ■ q) 2 ) , (C.5) 

where we used the well known expression 

W = -2 (8;5l - 815;) . (C.6) 

C.2 Calculation of the imaginary part of the reduced 
P — > l + l~ amplitude A 

To calculate the imaginary part of the integral A we use the rules of Cutkosky ( [69] ) where 
the off-shell propagators are replaced by 5-functions: 

\ -> -2iri5(k 2 ) 



k 

-> -2ni5{{q - k) 2 ) 



(q - ky 

Moreover we switch to the rest frame of the decaying particle P to rewrite its momentum 
as q = ( . The VMD part of course is equal to unity for on-shell photons. This leads 
to 
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We can now insert k 2 — k\ = k 2 and rewrite the (^-functions 

6(k 2 ) = Q(k )5(k 2 ) = ^-6(k 

S(( q -k) 2 ) = 5(m 2 P -2m P k ). 
After integrating over spherical coordinates this reads 



4 \p\ 2 \p 
which is exactly the desired result 



ImA = £-hx (- — § ) . (C.8) 



2/3 \l + p 

C.3 Derivation of the reduced amplitude -4.(0) 

We will give the derivation of the subtraction constant .4.(0) following the work of [13]. To 
evaluate the amplitude A(0) we first transform the integral from Minkowski to Euclidean 
space by k — > ik±. Now we will rewrite the VMD form factor in terms of a double Mellin 
transformation. According to e.g. [88] a Mellin transformation is defined as 

A(t) = M[A(s);t] = — [ A(s)r s dt (C.9) 
2m J 

with the inverse Mellin transformation given as 

POO 

A(s) = / dtA{t)t s ~ l dt. (CIO) 



o 



So, the double-mellin-transformed form factor VMD is given by 

vM Di e, (q - m = ^ dz * M (^) " . (an) 

Here A is the characteristic scale for the form factor, in our case the vector meson mass my, 
dz = dzidz 2 , the vector a = (<ti, (t 2 ) £ an d $(zi, z 2 ) is the inverse Mellin transform of 
the VMD form factor given as 

/*oo /*oo 

$(zi,22)= / dh dt 2 tf-Hf- l VMD{t u t 2 )- (C12) 
Jo Jo 

The integral f)3.5ip can be rewritten using Feynman parameters. Therefore the denominator 
reads (see e.g. equation (6.42) of [38]): 

1 



(k 2 ) zi+1 [(k - q) 2 } Z2+1 [(p- - k) 2 + m 2 ] 



3 \ „.«i„.«2 /-„.!-!> 



T(3 + Z 1 + Z 2 ) f -pr / ^ \ Q^Qg(Q^ , . 

y | j da, JM - ^ a, I — D]3+zi+22 (C.13) 



r(^ + i)r(i)r(* 2 + 1) 
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with D = aim 2 — a\a 2 q 2 . Applying this we can calculate the loop integral: 



2 f d A k (q ■ kf - q 2 k 2 



q 2 J tx 2 [k 2 + D] 3+Zl+Z2 

r( Zl + z 2 ) 



T(3 + Zl + z 2 )D 



ftq 



D 



-3 + 2+f m z - -q z ( Zl + z 2 



4 



(C.14) 



Using the Mellin transformed form factor and the loop integral (10. 14[) this leads to 



A(q 2 



(2m) 



dz 



-iR 2 



$(z 1 ,z 2 )(A 2 ) z ^T(z l + z 2 ) 

r(^ + i)r( 22 + i) 



X 



J U 3 i=1 daiS ( 1 - 



aim 2 — aia 2 q 2 ) zl+Z2 



-3 + 2 



2 (m 2 — \q 2 ) 



aim 2 — a\a 2 q 2 



(zi + z 2 ) 



(C.15) 



Because we are only interested in .4.(0) we set q 2 = 0. Using again equation (6.42) of 
to integrate over the Feynman parameters and setting £ 2 = m 2 /A 2 and z 12 = z\ + z 2 this 
reads 



.4(0) 



(2m) 



iR 2 



dz (e) 



-(zi+z 2 ) r(^i)r(z 2 )r(afi 2 )r(i - 2z 12 ) 



r(3 - z 12 ) 



X 



, -3 + 2z 12 )$(z 1 ,z 2 ) 

r( Zl )r( Zl ) 



(C.16) 



The integral can be solved, as in [13], with the help of the residues: 



A(0) = res Zr [Integrand#(0)]. (C.17) 
z r &n A 

The authors of [T3] followed the work of [5H] and found two contributions to the integral 
according to the residues 

a) z 2 + e = -ft, 

zi = -fr (C.18) 



b) z 2 + e = -ft, 

Zi = ft — (3. 



(C.19) 



C.3. DERIVATION OF THE REDUCED AMPLITUDE A(0) 



95 



They give the following contributions to the integral 

A(o) = A(o) + A(o), 



(C.20) 



A(o) 



^ (_!)«+? ^y+p+e T(-a -(3- e)r(l + 2(a + P + e)) 



a,6=0 



a!/?! 

x (_3-2(a + + e)) 



A(0) 



E 

a,6=0 



a!/?! ^ ' r(3 + /3) 



T(3 + a + e) 

Lr(-«)r(-/3). 

$(a — P + e, -a) 



(3-2/3) 



(C.21) 
(C.22) 



We need to give expressions for the Mellin transformed form factor. This can be done using 
Mellin transforms of derivatives as mentioned in Chapter 4.2 of [SB]- It can be proven that 
the following relation holds: 



M [x n F {n \x);s] 



-i)»^fw. 

MS 



(C.23) 



Here (x) denotes the derivatives of order of n with respect to the corresponding argu- 
ments of the form factor. Thus the following relation obviously holds: 

T(-a) 



dtr a - L F(t) 



dtt- a -H n F^ n \t) 



-1 



r(— a + n) 



(C.24) 



Setting n = a + 1 we can easily integrate over t: 

poo poo 

/ dtr a - l F(t) = (-l) Q+1 r(-a) / dtF^ a+1 \t) 
Jo Jo 



-l) a T{-a)F a {0). 



(C.25) 



In the last step we used F a (oo) = 0, which is obviously the case for our VMD form factor. 
If we extend this relations to two arguments of F, we can now give the following expressions 
for the Mellin-transformed VMD form factors 1131. 



T(-a)T(-P) 
*(z,-P) 

r(-a) 

Inserting (IC.26[) in (10.2 1 j) and ( 10. 22ft one gets the following result: 



■i) Q ^™^(o,o), 

POO 

-i) a / drViD (0 ' a) (i,o). 



(C.26) 



??! 



Mo) = E 

n=0 

MO) = E^t^T 



VMD n (0) , ~ n+e r(-n - e)r(l + 2(n + e)) 



r(3 + n + e) 



(3 + 2(n + e)), 



rfl + 2n) 



n=0 



nl 



r(3 + n)r(l -e + n) 



(3 + 2n) / dtt e VMD (n+1 \t). 



(C.27) 
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This can be expanded in e, and in the limit e — > we get 

r(l + 2n) 



= EH 2 )V 



n=0 



n)r(3 + n) 



yMD 2 (0) [2 + (3 + 2n) (ln4£ 2 - *(n + 3))] 
+ (3 + 2n)^ dtV^M.D n+1 (t) In t J. (C.2* 



To the lowest order of £ 2 one gets the desired result [13 



.4 (0) (0) = \ 



To the next order one finds: 



r POO 

31n£ 2 -- + 3/ dtVMD {1 \t)\nt 
2 Jo 



(C.29) 



-4 (1) (0) = 



VMD (1) (0) ^51n£ 2 + yj + 5jf cWML> {2) (t) In t . (C.30) 
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